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We study the energy minimization problem for an elastic interface in a random potential plus a quadratic 
well. As the position of the well is varied, the ground state undergoes jumps, called shocks or static avalanches. 
We introduce an efficient and systematic method to compute the statistics of avalanche sizes and manifold dis- 
placements. The tree-level calculation, i.e. mean-field limit, is obtained by solving a saddle-point equation. 
Graphically, it can be interpreted as a the sum of all tree graphs. The 1-loop corrections are computed using re- 
sults from the functional renormalization group. At the upper critical dimension the shock statistics is described 
by the Brownian Force model (BFM), the static version of the so-called ABBM model in the non-equilibrium 
context of depinning. This model can itself be treated exactly in any dimension and its shock statistics is that 
of a Levy process. Contact is made with classical results in probability theory on the Burgers equation with 
Brownian initial conditions. In particular we obtain a functional extension of an evolution equation introduced 
by Carraro and Duchon, which recursively constructs the tree diagrams in the field theory. 

PACS numbers: 05.40.a, 05.10.Cc 



I. INTRODUCTION: MODEL AND METHOD 

Complex systems, as well as systems with quenched disor- 
der, often respond non-smoothly, with jumps or avalanches, to 
a change in external parameters, as e.g. an applied field. This 
is seen as Barkhausen noise in magnets [JJ, earthquakes in the 
motion of tectonic plates |l2][3l, wetting of a disordered sub- 
strate f?,?!, dry friction p^l, cracks in brittle material dlSl, 
vortices in superconductors |9, 10| and many more. 

Quite generally, these systems can be modeled by an elas- 
tic interface pinned by disorder ll3l I9m8l |57^. In a previous 
work 1 19, 20 1 we have obtained the probability distribution 
of the sizes of static avalanches for an elastic interface in a 
random pinning potential. The interface is parameterized by 
a one-component displacement field u{x), where x is the d- 
dimensional internal coordinate. The interface is submitted 
to an additional quadratic well centered at w (e.g. a spring 
acting on it) and the ground state of the interface, u{x] w) ex- 
periences discontinuous jumps as the center-of-well position 
w is varied. These static avalanches, also called shocks be- 
cause of interesting connections with the Burgers equation in 
the limit d — 112114271 . are characterized by their size, i.e. 
S := u{x; w) where = J d'^x. In |fT9ll20l we obtained 
the distiibution P{S) from a combination of graphical and 
analytical methods, first at tree, i.e. mean-field level, valid for 
d > due, where the upper critical dimension is d^c = 4 for 
the usual short-range elasticity, and then to first order in a di- 
mensional expansion in e = 4 — d, by resumming all 1-loop 
corrections. This calculation was technically rather compli- 
cated as it required to sum an infinite set of diagrams, both 
at the tree and 1-loop level. Further their non-analytic depen- 
dence on w had to be extracted. Thus 09] |20l contains some 
amount of heuristics in extrapolating formulas from small mo- 
ments of P{S) to arbitrary ones, while the final structure is a 
relatively simple self-consistent equation. This suggests that 
a simpler method should exist. 

In this paper we present such a simple, complementary 
method. It is powerful and versatile enough to apply to many 
situations. It is extended in companion papers to (i) the depin- 



ning transition, for which the avalanche-size distribution was 
also studied numerically in |28 |: There we predict and mea- 
sure the distribution of velocities inside an avalanche |29-32|. 
(ii) Elastic objects where the displacements u{x) have more 
than one component ll33l |34| . This new method accounts 
for the (relatively) simple structures unveiled in ll20l . via a 
saddle-point equation and dressed propagators. It also allows 
to derive a more precise picture of the structure of avalanches 
around the upper critical dimension due- In particular we find 
that their statistics at due is given in the statics by a Brownian 
force model (BFM) which we study here, closely related to 
the so-called ABBM model ll35l that we also recently showed 
to describe interfaces at d^c near depinning ^291. 

In the second part of this article we make connection to the 
work by Carraro and Duchon |36|, as well as Bertoin f3T|. 
These authors use methods of probability theory to study the 
Burgers equation with Brownian initial conditions, which is 
the d — limit of the BFM for the interface. Their descrip- 
tion in terms of Levy processes is extended to interfaces and 
we unveil a new connection between evolution equations for 
these Levy processes in Burgers dynamics, and the mean-field 
theory for pinned interfaces. 

Our model is defined by the standard energy for a disor- 
dered elastic inteiface: 

n[u] =-HciM + [ Viu{x),x) (1) 
■^ciM = ^ / 9xx'Hx) - Wx]Hx') - w^'] , (2) 

^ Jxx' 

where g^x' = Jg 5^6*''^"^ ^ and the elastic energy kernel is, 
in the simplest case of short-range elasticity, = + m? 
(we denote = ^^y )- A quadratic external potential of 
curvature rn? and centered at Wx has been added and acts as 
a large-scale (infrared) cutoff. In all cases gq=Q = J^, gxx' = 
TO~^. V{u, x) is a centered Gaussian random potential with 
correlator 

V{u, x)V{u',x') ^ Roiu- u')5'^{x - x') . (3) 
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At finite temperature one considers the canonical partition 
sum Z = / 2?[w] e~^["l/-^ in a given disorder realization 
(sample). Disorder averages are denoted by and thermal 
ones by (•••). 

To study the statics of this model, one introduces replica 
Ua{x), a — 1, . . . ,n and considers the replicated action func- 
tional denoted iS/jo [u] = S[u] = S[{ua{x)}]: 

Z^= f n^^Kle-^I"! (4) 

a ab ^ 

The correlation functions of the disordered model are obtained 
from those of the replicated theory in the limit of n — > 0, 
implicit in all formula below. 

Let us now sketch the principle of the method, starting with 
a simple example. Since we are interested in probability dis- 
tributions of observables, we need to compute averages of the 
form 



ri^O \ / S 

where Ui{x) designates one of the n replica. One defines 



(5) 



(0M>5 



S\u] 



(6) 



and since for n = the denominator equals one, it can be 
dropped. Eq. (j5]l is a generating function from which one can, 
at least in principle, extract via Laplace inversion a probabil- 
ity distribution, here the distribution of the displacement field, 
i.e. Vlu] = (Yl^ 6{u{x) — u{x))) with a double average over 
sample and thermal realizations. Note that averages such as 
(|5]l, and their multi-point generalizations discussed below, are 
also frequently studied as generating functions of the distri- 
bution of the velocity field in turbulence. Burgers turbulence 
e.g. maps exactly to the present model at d = with time 
t = 1/to2 and velocity u itlTI l22l l25ti27l [36ti43l . 

We now recall a few basic facts from field theory. Let us 
consider T[u], the effective action functional associated to the 
action S[u]. Then the above average (j5| can be expressed as 



gJ a; ^ 



u^{x)-riu^] 



(7) 



where u^{x) extremizes the exponential, i.e. is solution of the 
equation 



(8) 



This property follows from the definition of the effective ac- 
tion as the Legendre transform of W[X] = (e-^^ 

from the relation W[X] + r[u] = J2a L Ua{x)X^. 

The calculation of r[u] can be performed in an e = due ^ d 
expansion around the upper critical dimension due- To lowest 



order in this expansion one replaces r[u] by the action S[u]. 
The corresponding calculation yields the tree-level result 



J^Ku\x)-S[u^] 



— Xx^al 



(9) 
(10) 



, A,trc 



which extremizes (|9|l. This precisely amounts to re- 
sum all tree diagrams in the perturbation expansion in the 
non-linear part of the action Shq, i e. in the disorder i?o, also 
known as the mean-field calculation. This is discussed in sec- 
tion [III] and appendix |D] For the problem at hand it gives the 
correct result for probability distributions for d > due, if the 
renormalized disorder R is used in the action, rather than the 
bare one Rq, as discussed in |!20| and again below. The corre- 
sponding action Sr is called the improved action. The precise 



definition of R is recalled in Section IV A 2 and useful equiv- 
alent definitions can be found in Sections II and III of Ref. 
II20I (with which present definitions and notations aim to be 
consistent). 

Here we start with the tree calculation in section |Il] by first 
defining the proper observable to compute, a generalization of 
(|9]l, and deriving the saddle-point equation which resums all 



tree diagrams. In the following section III we give a graphical 
derivation and illustration of the saddle-point equations. We 
then compute in section IV r[u] to first order in e, and analyze 



the resulting saddle-point equation, from which the avalanche- 
size distribution to 1-loop order is obtained. In Section [V] we 
study a simpler model where the force landscape is a Brow- 
nian motion, and we make the connection to Levy processes 
and Burgers equation. In Section |VI] we derive the general- 
ized Carraro-Duchon equation which encodes the mean-field 
theory of interfaces. 

The appendices contain details and extensions: In appendix 
[a] we study a non-uniform deformation Wx ■ Appendix [B] 
derives useful formulas for the diagonalization of replica- 
matrices. Appendix[c|calculates Fi [u, v]. Appendix[D]gives a 
diagrammatic interpretation of the loop corrections. Appendix 
[E] contains a detailed derivation of many-point correlations in 
the BFM. Appendices |F| and |H] recall the derivation of the 
Carraro-Duchon formula and its connection to the exact RG 
equations. In Appendix |l] we discuss the (near absence of) 
loop corrections in the BFM model and we prove that it is an 
attractive fixed point of the RG. Finally in Appendix |j] we re- 
call how the statistics of shocks depends on their correlations. 



II. TREE-LEVEL (MEAN-FIELD) CALCULATION 
A. Avalanche Observables 



Let us recall the avalanche observables introduced in 
Unless stated otherwise the considerations below are valid in 
all generality (i.e. beyond mean field). 

At T = 0, the minimal energy configuration of the interface 
u{x; w), and its center of mass u{w) := L~'^ u{x; w), ad- 
vances by jumps as the position of the center of the parabola. 
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Wx (taken uniform for now), is increased: 

u{x; w) — Si{x)9{w — Wi) 

i 



where 



(11) 
(12) 



i?"'(0+) 



(18) 



Here Wi is the position, and Si := J^Si{x) the total 
size of the i-th shock. One defines p{S) :— pa P{S) = 
J2i ~ Si)S{w — Wi), where p{S) is the shock-size den- 
sity, P{S) the shock-size probability distribution (normalized 
to unity) and po the total shock density. From u{w) = u> it fol- 
lows that Po (S) = L"^ whenever all motion occurs in shocks 
(which is the case here in the limit of interest, m — ^ 0). We 
denote size moments as {S'f) :— J dS SPP{S). It was also 
shown in H20J that for 5*0 ^ S, P{S) takes the general form 



is the characteristic avalanche size introduced in Eq. (14 1. 
Taylor-expanding Eq. (17 1 in A, Z{X) = A + S'^A^ + . . . , 
and comparing to the definition ([15]), allows to identify Sm = 
i(52)/ (5), also stated in Eq. ([^ 
Inverse Laplace transforming t 



17 1, one obtains the 



mean-field avalanche-size distribution , valid for d> 



Ptr 



= (s) = Pmf(s) = 



1 



-s/4 



(19) 



PiS) 



(S) 



p{S/Sm) 



(S^) _ i?'"(0+) 
2(5') ^ m4 



(13) 
(14) 



Here S„i ^ m^'^^'» is the scale of the large avalanches, and 
So <C Sm a microscopic cutoff. The function p{s) is universal 
with ds sp{s) = 1 and ds s^p{s) — 2. 

As detailed in [20 1 the shock-size moments can be extracted 
from the generating function: 



Z{X) := L-'^a„c^i''[«(«'/2)-"'-"(-«'/2)] 

_ (e^^) - 1 - X{S) 
{S) 



w=0+ 



(15) 



We now recover these results, and more, by introducing a 
method which does not use a graphical expansion. 



B. Saddle-Point Equation 

Here we show how to evaluate, at tree level, the slightly 
more general generating function for the joint probability of 
the field u{x; w) at two "points" and —Wx/'i. It corre- 
sponds to moving the center of the parabola from —Wx/^ to 
Wx/2. While this is not the most general non-uniform move, 
its symmetry simplifies the analysis below. For future con- 
venience, we denote here the full disorder correlator (which 
contains loop corrections to all orders), by R, i.e. we consider 
the improved action Sr. This is an improved tree approxima- 
tion, i.e. it is the sum of all tree diagrams in R; in Rq it is 
the sum of all tree diagrams plus those loop diagrams in Ro 
correcting R itself. 

Generalizing Eq. (j9]) requires to introduce two sets of n 
replicated fields denoted Ua,Va, a = 1, . . . ,n, subject to the 
same disorder. We find that^ 



This formula follows from the fact that in a small window of 
width w > the probability that there is a shock is pow, in 
which case the field u{w) — w jumps by S. We used that 
due to statistical translation invariance, we can without loss of 

generality consider the interval ] - w/2, w/2[. Equation ^ where (dropping the superscript A): 



p/_^ Xj;[u(x;w/2)-w^-u(x;-w/2)] 



(20) 



can be compared to (jSj) with a uniform A^ = A: however, 
while (jsj) encodes only the one point probability of u{w) (say 
at u; = 0), Eq. ( [T5] l depends on the two-point joint probabil- 
ity distribution of the field u{w) at two values of w (denoted 
—w/2 and w/2), as required to study shocks. Furthermore, 
in d = 0, 'm?{w — u{w)) is the velocity field of the decaying 
Burgers equation at space point w 121112211251 and Eq. ( 15 i is 



Sx[u,v\= I \x[ui{x) - Wx ~ vi{x)] 

J X 



Wx 

2 



■ . ,x Wx' 
Ua(x ) - — 



thus the generating function of the distribution of velocity dif- 
ferences at two points in space distant by w. We retain below 
this terminology of p-point distributions. 

We now recall the results obtained at mean-field level in 
II20I by resumming the tree diagrams. There it was shown 
that Z{X) := A + Z{X) satisfies the self-consistent equation 



Vaix) + ) + I 

J [R{ua{x) - Uh{x)) + R{Va{x) - Vb{x)) 



+ 2R{ua{x) ~ Vb{x)) 



(21) 



Z(A) = A+^^Z(Af . 



(16) 



This quadratic equation is easily solved, 

Ztree(A) = Zmf(A) = ^ (l " - 45„,a) , (17) 



' These formula correspond to an infinite density of avalanches. For discrete 
displacements u, as illustrated in Appendix|j] po is finite, and Z{X) is cut 
at large negative A, equivalent to a small S = Smin for P{S). This is 
further discussed in (20|. 

^ Here u{x; w) is a functional of the field . For simplicity however we do 
not use the square bracket notation. 
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Ua{x) and Va{x) are extrema of Sx, i.e. solution of 



with 



TX 



= gj, [uaix') - ^ (22) 
- ^ 51 - uc(x)) + R'{ua{x) - v,{x))] 



together with a similar equation for v, obtained by 
{u, A, w) -H> (v, —A, — ui). One can also write the functional 
derivative of (|20ll. 



(23) 



2r 



Ua{y')-Wy'-Va{y') 



Due to the saddle-point equations p2] ), only the explicit de- 
pendence on w appears. 

By parity, the solution of the saddle-point equation satisfies 
Va{x) = —Ua{x), which we use from now on. Since only 
replica 1 is singled out, we look for a (replica symmetric) so- 
lution where all replica a 7^ 1 assume the same value. We 
denote 



Uaix) = Ui{x) -TU{x) , 07^1, 



(24) 



which at this stage is just a definition. The saddle-point equa- 
tions for the two functions ui{x) and U (x) become 



1 



(25) 



T 



R'{TU{x)) + R'{2ui{x) - TU{x)) - B! {2ui{x)) 



KT^T I ^ g-^,U{x') - -R'{2u,{x) - 2TU{x)) (26) 
[2R'{2ui{x)-TU{x)) - R'{2ui{x))] . 



Note that the second equation has been obtained by subtract- 



ing in (22 1 the equation for a ^ 1 from the one for a = 1. 
In both equations the sum over replica indices has been per- 
formed using (|24|, i.e J^c ^("c) = F{ui) + {n - l)F{ui 



TUi), and then setting n — )• 0. We have also used that R'{u) 
is an odd function with i?'(0) = 0. Once these equations are 
solved, the solution can be used to compute the generating 
functions: 



e 

■tree 



(27) 



-Sx:=-Sx[u,-u] (28) 
= / Xx{2ui{x) - Wx) 

U{x) {2ui{x') - w^^ + TU{x)U{x') 

2R{0) - 2R{TU{x)) + R{2ui{x)-2TU{x)) 

+ R{2ui{x)) - 2R{2ui{x)~TU{x)) 



+ 



+ 



9xx' 



J^2 



The limit n = has been taken everywhere. This result is 
equivalent to the graphical summation of all tree diagrams, 
in terms of either i?o or R, depending on whether bare or 
renormalized perturbation theory is used. It is valid for any T 
and Wx, hence in principle it allows to compute at tree level a 
rather general 2-point correlation function of the field u{x; w) 
at any temperature T. 

In the absence of disorder, the saddle point is ui(x) — 
Wx/2 — TU{x) — T J^,gxx'Xx', and one obtains —S\ = 
T ^1 gxx'XxXx' ■ The tree formula d27b is then exact and 

elatei 



corresponds to two copies with uncorrelated thermal fluctua- 
tions. In presence of disorder, but for = 0, one must have 
Sx = 0. This is indeed the case, as the saddle point is then 
U{x) — and ui{x) — ^ ~ 0. When there could be several 
solutions to the saddle-point equation, the correct one should 
reduce to that one in the small-A limit. The saddle-point so- 
lution can also be obtained order by order in A from perturba- 
tion theory, i.e. a well-defined expansion of ui{x) and U{x) 
in powers of A must exist. 



C. r = limit of tlie saddle-point equations 

We can now study the system at T = 0. Then R{u) is non- 
analytic, more precisely it exhibits a linear cusp in its second 
derivative, _R"'(0+) > 0. This cusp is related to the second 
moment of avalanche sizes, as shown in |20|, via i?"'(0+) = 
m*^ (S^) /2 (S). We will recover this relation here. 

In the T — ^ limit we obtain a consistent solution assuming 
that U and mi are going to a finite limit, as we show now. 



Expanding p5] l and ( 26 1 in powers of T we obtain to lowest 
order: 



ui{x') 



9 XX 

= 

.— 1 TT( TDIII 



- [i?"(0)-i?"(2wi(x))] U{x) 

(29) 

gxx'U{x') - R"'{2ui{x))U{xf = A, (30) 



The generating functions are, omitting the thermal averages 
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. . ), since we are studying T = 0, 

\j.lu{x-.w/2)~w^~u{x-,~w/2)] _ Q-Sx 
Jy' 



-Sx = 



ui{x) - 



2 A. - / g-,\uix') 



(31) 



(32) 



(33) 



where (29 1 has been used to simpHfy 5a. These formula are 
valid for arbitrary Wx ■ 

We now analyze these equations in several cases. 



D. Uniform case: avalanches of center of mass 

Let us start with the simplest case of both = X and 
w.r = w uniform, corresponding to the generating function 



( 15 I. Then ui{x) = ui and U{x) — U satisfy 



m^{ui - w/2) + [R"{0) - i?"(2Mi)] U ^0 

m^U - E!" {2ui)U'^ = X 

L~'^dyj e^-t''["(™/2)-'i'-«(-™/2)l'"'' 

= [-X + m'^U]e-^'' 
-5a = L'^{2X-m^U) (ui - ■ 



(34) 
(35) 

(36) 
(37) 



These equations can be studied for any w. 

We now consider the limit w from which avalanche 

observables can be extracted. We look for a solution of the 
form 



which implies that S\ = 0{w). Hence 

-tree 



w=0+ 



Assuming y > we obtain from p4\ and (35 1 



m^U - B!" {Q+)U^ = X 



(38) 



-X + rri^U . 
(39) 

(40) 
(41) 



Comparing the definition ( [T5| and the result ( 39 1 we see that 
we can identify 



Z{\) = A + Z(A) = m'^U 



(42) 



Our self-consistent equatio n (|4T] ) is indeed the same as the one 
obtained in |20|, namely (|16|l. Its physical solution, which 
vanishes at A = 0, is m?U — Z{X) — Zmf(A) given in ( 17 1. 
Note that 



1 



y 



(43) 



is indeed positive for this solution. The breakdown at A > 
1/(4S'.,„) signals that the Laplace transform of P{S) does not 
exist beyond that value of A, due to the exponential tail at large 
5 in 



E. Non-uniform case: Local structure of avalanches 

To obtain spatial information about avalanches one may 
consider both Aa; and non-uniform. We specify w^. = 
wf{x) and vary w from w — to w = +^ for a fixed 
f{x). In the limit of small positive w we look for a solution 
of the form 



w 



= y{x)^ + 0{uP) 
Again one finds 5a = 0(w), hence 

Q A^[u^(to/2)--!i>^-ji^(-u)/2)] 



tu=0+ 



(44) 



(45) 



/ 9-x,U{x') 



The field U{x) satisfies Eq. (30i, i.e. 

g-^,U{x')~R"'{y{x)w)U{xf^Xx 



(46) 



Here we assume that y{x) > 0, a more general discussion is 
given in Appendix [a| Then in the limit of w 0+, one can 
replace R'" {y{x)w)-^ i?"'(0+) and (46l becomes a closed 
equation for U (x), independent of y{x). 



9~^,U{x')-R"\Q+)U{xf 



(47) 



This is a classical equation for a cubic field theory, which ad- 
mits instanton solutions, from which local size distributions 
can be extracted, as discussed in |20|. Here we will not study 
again these applications, but simply make contact with the no- 
tations used there. For that purpose we identify 



Zx{\) = / 9-x'U{x') , (48) 

J x' 

in terms of which the self-consistent equation becomes 



^.(A) = A, + i?"'(0+) 



QxXx 9xX2 ^Xi 



(49) 

These are Eqs. (204) (in unrescaled form) and (F8) of |20|. 
The space-dependent generating function can then be written 

as 



— ^a(A) — Ay := Zy{X) . 



(50) 



Note that a rescaled version of U{x) was denoted Y{x) in 
I20J. 
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Zy{X) is connected to local avalanche-size distributions. 
Assuming that for fixed L as w 0^ the probability of a 
shock during a change of Wx is plw, we can write from Eq. 



Pb 



f UL^.s^ 



1 



(51) 



The subscript / reminds that we use a non trivial f{x). In 
addition, one can show that Ux{w/2) — Ux{—w/2) — 
hence Pq {Sx).f — 1, which also implies Pq (S) ^ — L'^. Note 
that the size distribution, whose Laplace transform is given by 
(51 1, a priori depends on the function f{x); the case f{x) = 1 



was studied in 1201 . 



F. Multi-point correlations of center-of-mass displacement 

1. Discrete version 

We now indicate how to compute, at tree level (i.e. as a sum 
of all tree diagrams), the correlations of the center-of-mass 
displacement field u{wi) — wi at an arbitrary number of dis- 
crete points Wi. For this we introduce a generating function, 
parameterized by A^: 



gi^'E, ^.["("-O-w-il (52) 
The fields Uai are solutions of the saddle-point equations 

m,'^{Uai - m) - Tf;^ R'{Uai - Ucj) = TXiSal ■ (53) 



As above, we look for a replica-symmetric saddle point Uai — 
Ui for a 7^ 1. Define Ui := uu — TUi, and subtract the 
equation for a = 2 from the one for a = 1: 



m'^{uu - Wi) 
1 



(54) 

= TA,; 



j 



+R'{uu - u^o -TU,)- R'{uu - u^j - TU, + TU^) 
= TA, (55) 

Note that in the second equation the terms i — j can be ex- 
cluded. Taking the limit T — > 0, we find 



■{uii -w^) + Y^ R"{uu - uij)Uj = 
j 

m'^U, + ^ R"\uu - ui.j)U,U, = A,; 



(56) 



One must solve these equations for Ui and uu and insert the 
result into 

ei^^ET^TR^'JO^*'''"' = e^''E.(A.-i™=c/,)(«i.-™,) ^ (57) 

which has been simplified using the saddle-point equations. 
Note that one recovers the 2-point equations ( [35| from the 
solution U12 = —uii and U2 = —Ui valid for A2 = — Ai. 



2. Continuous version 

It is instructive to perform the same calculation in a contin- 
uum framework. We again restrict to the center of mass and 
compute 



-tree 



(58) 



for a test function \{w), depending on w but uniform in space. 
Here and below f ~ dw. We now introduce^* replica 
fields Uai Ua{w) and then extremize the action: 

L-'^Sx = [ X{w)[ui{w) - w] 



+ ^ / Y.^{ua{w)-u,{w')) . (59) 



We will not repeat all the above manipulations. The saddle- 
point equations w.r.t. Wa(w) in the limit T = lead to 



jn^luiiw) -w]+ / R"{ui{w)-ui{w'))U{w')=0 

J w' 

m^U{w)+ I R"'{ui{w) - ui{w'))U{w)U{w') = X{w) 



(60) 



Its solution is inserted into 



X{w) ^ —U{w^[ui{w) -w] . (61) 



The general analysis of these multi-point correlations, and 



their discrete analog ( 56 1 is left for the future. Below we study 



them, in Section |V] and Appendix |E] for a simpler model, 
where R"'{u) is a constant, and in Section VI G for periodic 
disorder We also discuss, in Section VI another powerful 
method to generate multi-point correlations. 



' The generalization with space dependence would involve A^; (w) and 
replica fields Ua (x, w). 
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G. Displacement correlations for finite w: Absence of 
correlations for d > due 

The two-point correlation function of the center-of-mass 
displacement. 



qL'^X[u{w/2)~-w-u(-w/2)] _ g-5A (-g2) 

can also be evaluated within tree level at finite w > hy 
setting Ml = yw/2, and solving the saddle-point equations 
( [35] l for y and U, denoting y„, and 11^ these solutions. This 
can be done, e.g. order by order in w. Here we give the small- 
w expansion up to 0{w'^). It is convenient to introduce the 
rescaled (dimensionless) variables w and A, 



W = m SrnW , A — \/Sm ■ 

Defining = m^SmUw, one finds 



(63) 



{mL'f 

Z = ^trcc(-^) = 



2(1 - 2Z) 



i?""(0) 



2^(1-Z) i?"'(0)i?(5)(0) 



(1-2Z)5 6(1-2Z)3 



1 -4A 



(64) 
(65) 



More generally, one can introduce the dimensionless rescaled 
renormalized disorder correlator. 



R"{u) ^ Adm'~^'^R"{um'^) , 



(66) 



where Ad — l/(e/2), and for short-ranged disorder e = 4 — d 
and 



1 



(67) 



is an amplitude; for details see ifTT! . and for generalization to 
long-range disorder lfT6l l44l l45l . It is known that as m — > 0, 
R"{u) goes to a fixed-point function in any d, measured in 
B6l . Then, in any d, and within the tree approximation, 

" = " ^^(^ " ^^^^ 

takes a scaling form, obtained by eliminating y = y^ and 
Z = Zyj in the rescaled saddle-point equations: 

(, - i)AdR"'io)^ + R'm-R'W-jO)^)^ ^ ^ 

^ ' ^ '2 i?"'(0) 



z - 



R"'{yAdR"'{Q)w) ~ 



R"'{0) 



Z^ = \ 



(69) 



While these equations for the tree approximation can be writ- 
ten for any d < d^^, they are expected to become exact as 



due and in d > due- For d 



e, e > 0, the rescaled 



disorder R"{u) is uniformly of order e. In (69 1 the argument 



of the functions is 0(e) which justifies an expansion of the 
system for small w. This is because the scaling variable is 

w = AdR!"{Q)m-'^w = Smm'^w . (70) 

Hence near d = due one can focus on (|64]). Since 

m-'i?""(0) = Arfi?""(0) = 0(e) (71) 
m-2^i?"'(0)i?(5)(0) - O(e') , (72) 

we arrive (in terms of the unrescaled variables) for w > at 
the result 

gL''A[«(tu/2)— !i)-ti(-ii;/2)]''°° _ (,L''wZ{\)+0{e) ^y^) 

The interpretation of this result is that, for d > due the in- 
crements in the displacements in the center of mass become 
uncorrelated. It will further be discussed in Section IV Bl in 
terms of Levy processes. 

Below due, in d = due — we expect correlations. They 
only exist on a distance w — 0{1), i.e. w — 0{e)m~'', a 
very small layer as e — > 0. The above result (j64| allows to 
compute the first correction in e; however one may also get 
contributions at one loop of the same order This calculation 
is performed elsewhere. 

III. GRAPHICAL INTERPRETATION 

Here, we sketch a short graphical interpretation of the 
mean-field saddle-point equations. We work in d = for 
simplicity. The results apply to the center-of-mass variable in 
any d, after restoring the necessary factors of L"^. For easier 



comparison with section VI we use the notations 



A(u) = -R"{u) 



t = 



1 



We define 



„A[«(to)-M(0)] 



(74) 



(75) 



Thus Z(A, w) = Z(A, w) — wX is the generating function of 
the connected moments [u{w) — u(0) — w]" — fif^"-' (w), 
which in Eq. (41) of |20| were called the Kolmogorov- 
cumulants. In EOl . a graphical derivation of the recursion 
relation for the 0{w) term Z{X) was given, noting 



Z{X,w) ^ Z{X)w + 0{w^ 



(76) 



At tree level, Z^^{X, w) is the sum of all connected tree dia- 
grams. The generating function Z^^{X) as a function of the 
bare action can be written as a sum of particularly simple tree 
graphs, namely the ones of the form 



ZMF(A)=^^ 



i 



i 



(77) 
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where dotted lines represent the bare disorder Aq. Each graph 
represents correlations of the X[u{wi) — Wi] fields at different 
points Wi (external lines on the top, each coming with a factor 
of A). The linear operator IC identifies Wi with w or 0, in order 
to build the Kolmogorov cumulants. The disorder vertex on 
the bottom contributes a factor of Ao(w) — Ao(0) to Z(A, w), 
which has been expanded to first order in w to obtain Z{X), 
hence it must be counted as Ag(0+) in (77i. Equivalently, 
one can include a factor of duj\w=o in IC. For more details of 
these graphical rules see f2n\, section V.C. They are also fur- 
ther used below in SectionlVljand appendix|D] The improved 
generating function Z^^- '^{X) is the same sum of tree dia- 
grams with A at each vertex, hence, if reexpressed in terms of 
the bare disorder Aq, it is now a sum of graphs of the form: 



(78) 

It contains all loop corrections to the 2-point disorder vertex, 
while loop corrections to higher-point disorder vertices are ne- 
glected. Explicit formulas for low-order contributions can be 
found in EUl . 

The mean-field self-consistency equation for Z{X) = 
ZMF.^B.^X) reads 



Z{X) = X-t'^A'{0+)Z{Xy 



(79) 





It is graphically written as ||20] 
Z{X) ' ^ 



As indicated, the blob denotes Z{X) itself, while the low- 
est disorder vertex counts as a A'(0"'") and the lines enter- 
ing the two blobs from below do not come with differentia- 
tions'*. This self-consistent equation yields the desired sum of 
tree diagrams with only one lower disorder vertex, i.e. a term 
A(w) — A(0) expanded into wA'(0+), which is sufficient to 
obtain the 0{w) part in (76 1 as explained in 



We now want to construct a recursion relation for Z(A, w), 
which yields its complete w dependence. We thus need to sum 
all tree diagrams. To generate them, it seems natural to write 



Z(A,w) 




— Xw — 




(81) 

Now the lower vertex is A{w) — A(0) and the line entering a 
blob Z(A, w) from below acts as derivative w.rt. w. 



* Note that in [201 a factor of | A'(0+)| = - A'(0+) has been absorbed in 
A, whereas here it is explicitly written, resulting in a seemingly opposite 
sign. 



However, a new difficulty arises: One may have two or 
more lower vertices, as e.g. in 



(82) 



Unfortunately, the self-consistency equation ( 8 1 1 is then in 



correct, as it leads to an over-counting since there are several 
ways to construct the same graph. Fortunately this can be 
corrected. Let us explain the source of the problem, and its 
correction on an example: 



J i 1 



1 i \ 1 



(83) 



On the l.h.s. we have plotted the contribution to Z(A, w), with 
the correct combinatorial factor The first two terms on the 



r.h.s. appear in a recursion relation of the form (81 1, plotting 
the added vertex in Eq. ( [ST] ) red (fat in black-and-white). This 
leads to an over-counting, which can be corrected by subtract- 
ing the last term, which has two marked (red/fat) vertices. 
In the case of three lower vertices, the recursion reads 



+ 



(84) 



Now pairs of lower vertices are subtracted, leading to a can- 
cellation of all terms, and consequently the triplet of lower 
vertices has to be added. This can be generalized, replacing 



the self-consistency condition (81 1 by 



(85) 



The dots around the first exponential function indicate that 
the derivatives act only on Z(A, w), not on A{w). This can be 
written as 



Xw = Umiyln (e-[A(-)-A(o)]t^O„a„giz(A,n)^ 



(86) 



The limit of selects the tree diagrams, and the In selects 
a single connected component. Note that e''['^("')~^('^)l* 
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is defined by its series expansion in t, thus the limit is 
done term by term in the expansion in powers of t. 

To proceed, we observe that independent of i' and for all 
functions f{u) which are infinitely differentiable^ 



(87) 



^47riy[A(w) - A(0)]t2 
Inserting this relation into ([86| yields 



due 4.[A(™)-A(o)it^ 



Xw — lim In 



^47ri/[A(w) - A(0)]i2 



(88) 



In the limit of — to be taken here, the integral is domi- 
nated by its saddle point, and we get 

wX ^ Z(A, u) + [A(0) - A(w)] f [9„Z(A, u)]^ (89) 
u = w + [A{w) - A(0)] a„Z(A, u) . (90) 

The new variable u has to be eliminated between these two 
equations. Note that it is an independent variable not to be 
confused with u{w). 



We remark that passing from ( [86| to the self-consistent 
set of equations 



pQf is a quite common feature in tree- 
resummation problems. It also appears in the large- resum- 
mation for the disorder itself, where the links are the 1-loop 
momentum integral, and which therefore are termed cactus- 
diagrams, see [47 , 48 1. 

Below we devise another method to compute Z(A, w) and 
we have checked to high orders, ^ A^^'', that (89 1 reproduces 



the solution Eq. (210i, e.g. (212 1 as well as the lowest-order 
Kolmogorov cumulants (62)-(66) in ||20]| . 

One needs the derivative obtained by deriving (90i 
w.r.t. w. 



du 
dw 



l + 2t2A'(w)(9„Z(A,it) 
l + 2t2 [A(0)- A{w)]dlZ{X,u) 



(91) 



Deriving Eq. ( 89 1 w.rt. w, and using (j9T|i yields the astonish- 
ingly rather similar equation. 



A = a„Z(A, u) + A'{w)t^ [5„Z(A, u)f 



(92) 



We now make contact to the results given in Eq. ( [34] l. Since 
u and w are simply variables, and we are more interested in 
Z(A, w) than in Z(A, u), one can exchange their names, to 
obtain a second set of equations. 



uX = Z(A, w) + [A(0) - A{u)] f [a^Z(A, w)] 
w~u^2t^ [A{u) - A(0)] 5^Z(A, w) 

X = d^Z{X, w) + A'{u)t^ [a„Z(A, w)f . 



(93) 
(94) 
(95) 



^ The formula is valid for v{A.{vj) — A(0)] > 0, i.e. the integral should be 
evaluated for f < 0, and the limit to be taken is ^ — ^ 0~ . 



Eq. ( 95 I is redundant, or can be used instead of ( 93 1. We have 
explicitly checked that up to order t^, both expressions are 
correct. 



Graphically, the interpretation of Eqs. (89h92i is rather 



different from that of Eqs. (93 1-(95 i. To see this, recursively 
replace u in Eq. (|89]l by its value given by Eq. ( [90| . This 
yields a perturbative expansion in t of Xw, which can be read 
as a self-consistent equation for Z(A, w). Graphically, it con- 
tains links made out of [A(w) — A(0)] t^, which end in ver- 
tices made out of Z(A, w). If n links enter into such a vertex 
Z(A, w), it means to take n derivatives. 

The picture is different when replacing recursively in Eq. 
( 93 1 u by its value given by Eq. (j94]|, thus constructing again 
a perturbative expansion in t for Z(A, w). Apart from a single 
term Z(A,u'), all other terms are proportional to powers of 
9u,Z(A, w), and no higher derivative of Z(A, w) appears. The 
objects with more derivatives are [n — 2)-nd derivatives of the 
disorder [A(w) — A(0)]t2, which have n outgoing lines which 
end either in a disorder, a 9^Z(A, w) or A. This will in more 
detail be discussed in |33 1. 

Identifying 

a„Z(A,?i;) = m^C/ (96) 
u = 2ui (97) 



shows that Eq. (95 i is equivalent to Eq. (35i, and (94i to 
(j34]). We still have to check expression We know that 
Z(A, It;) Z(A,?i;) - Aw = ~Sx. Multiplying Eq. (94 1 with 
A, and adding the result to ([93]l gives 



Z(A, w)-Xw^ [A{u) - A(0)] [d^1{X, w)f 

-2t^X[A{u)- A{Q)]d^Z{X,w) (98) 
The r.h.s. is nothing but —S\, given in Eq. ( [37] i. 



IV. BEYOND MEAN FIELD: 1-LOOP CALCULATION 
A. Simpler example: 1-point probability 

1. Perturbation around mean field 

We Start with the simpler case of a 1-point probability, e.g. 
as given by Eqs. (|7| and ([8]). As we explain in detail below, 
in the dimensional expansion around mean field, the effective 
action can be written as 



r[ii] = 5M + riM 



(99) 



Here S[u] = SrIv] is the improved action, and ri[7i] is 
"small" in a sense to be specified below. Hence we can expect 
that and u^'^'^°°, the solutions of (jsjl and ( lOi, are close to 
each others. Schematically we write 



^/r A.treen 



A 
A. 



(100) 
(101) 

Inserting into Eqs. (j7]i and (jHll and expanding to lowest order, 
i.e. in the differences — ^ 0{Ti) and Fi, we find 



(102) 
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Hence to compute this generating function to lowest order 
around the tree result we only need to evaluate Fi at the tree 
saddle point. 



2. Effective action for the pinned interface 

The (replica) effective action V[u] associated to the (bare) 
action Q of the pinned interface takes the form 

rM=5eiM-5] (103) 

p—2 ai ,. . . ,ap 

Here Sc\[u] — ^Tlia-^c\[uo\ arises from the elastic and 
quadratic well energy, defined in (|2|. The p-replica terms, 
S'^P\ are the p-th cumulants of the renormalized disorder*. 
The local part (i.e. Ua{x) — Ua) of the second cumulant, 
5*^^^ [wa, Uf,] = L'^R(uab), Uab — Ua — defines the renor- 
malized disorder correlator R{u). 

Let us now consider T = 0. As m -H- 0, R{u) flows to a 
fixed-point function R — 0{e) where e — d^c — d. In general 
r[u] can be computed in an expansion in powers of R, i.e. of 
e. Forp = 2 one has ifTTllBl : 



"'^ ab 



(104) 



1 



9 , (gly - 5,ygl)R"{uabix))R"{uabiy))+0{R') , 

^ Jxy 

where here and below we denote R"{u) := R"{u) — R"{0). 
Purely nonlocal parts are thus O(e^) and higher. For p > 3, 
each S'(P) is of order 0{RP) = 0{eP). They were computed 
previously ll20ll49l : the result can be summarized by Eq. (99 1 
with 



1 



1. 



T,[u] = -Trln {g-'6ab - W^,) - -Trln (g'^S^b - 



ab) 



(105) 



W:b^,y = ^S^yiSab R"{^ac{x)) - kR" {uab{x))] (106) 



where /„ — Xc ffl' ■ ' 105 i has the usual expression of a 1-loop 



effective action, ^Tr In 5", apart from a subtraction of the 1- 
loop graphs leading to p + 1 replica terms, proportional to T 
(second term on first line) and of the p ~ 2 part, akeady taken 
into account in (104i, thus expressing ri[u] as a functional 



of the renormalized instead of the bare disorder. (Note that 
as written Fi also contains the bilocal 0{e^) part of p = 2 



in Eq. ( 104 1). Upon expanding the Tr In, which acts both on 
replica and space indices, the S^-^^ are recovered to OlW^); 
see e.g. formula (113) in 120 J . Note that the two 0{W) terms 
cancel (using n = 0). 



3. 1-loop probability-distribution of displacements 

Let us compute the 1 -point generating function to lowest 
order beyond mean field. 



„-ri[{u„(x)}] 



(107) 



Here Ua (x) satisfies the tree-level saddle-point equation 

J x' ^ ^ 

(108) 

Supposing that there are exactly two different fields ui and Ua 
for a > 1, this gives 

9xx'\Mx') - w^.] + ^R'{ui2{x)) = TK (109) 



9xx'[<^2{x') - W,,] - -R'{u2l{x)) = (110) 



Note that in the first line we have dropped the term ^ R! (0) = 
0, and the sign change comes from the factor of 1) replica. 
On the other hand, from the sum in the second line, only the 
term c = 1 survives, while all other terms are ~ i?'(0) = 0. 
As in equation p4|, we now look for a solution Ua — 



Eqs. ( 109 1 and \\\Qi) become 



ui — TU for a 7^ 1. Eq. (109i and the difference between 



TXx = fg:^,[ui{x')-w,,] + ^R'{TU{x)) (111) 



Ux ^ gxx'K' (112) 

J x' 

In the limit of T — > 0, the first equation is expanded as 

ui{x)=w^^R"iO) [ g,^.U{x') + 0{T) 



R (0) / / gxygyA. + OiT) . (113) 

J y J z 



This gives the tree contribution at T = 0, 

-tree 



T=0 



(114) 

It is a Gaussian distribution for the displacements at tree level 
(recall i?"(0) < 0). 

Let us now compute the 1-loop corrections to mean field. 
Here we restrict to uniform = X and uniform Wx — w. 
Since we consider a 1 -point function, w = can be chosen. 
The saddle point is uniform Ua{x) = Ua, and 



' Due to statistical tilt symmetry, the term p = 1 is a constant, dropped here. 



^^trln (g^^Sab + M^,) 



tr In {g-^ Sab + M^,) 



+ ^tr[(Mif-(MO)^] 

^[Sab [R"iUac)-R"{Ual) 



(115) 

KR"{Uab)}in6) 
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Here tr refers to replica indices, and we have used the saddle- 
point properties denoting by c any index with c 7^ 1. More 
explicitly, is a replica matrix, where replica 1 is singled 
out, and which is symmetric in the other n — 1 replica a 7^ 1. 
It thus has (for any n) four distinct components (denoting with 
a, b indices different from 1): 





^^[R"iTU) + in-l)R"m 


(117) 






(118) 






(119) 




^^[i?"(o)~i?"(rc/)] 


(120) 




= ^^i?"(0) 


(121) 



The field U isU = Xjm}, see Eq. ( 1 12 1. The diagonalization 
of this matrix is performed in Appendix [B] The eigenvalues 
are given in Eqs. ( |B5| l and ( |B6| l. In the limit n they are: 

(i) a (—2) -dimensional space with eigenvalue [i — Mc, 
since A/^ is independent of k, its contribution cancels be- 
tween the two first lines of (11 ISb. 



(ii) a 2-dimensional space with eigenvalues given in (|B6| l. 
Since /2'" = Af{\ +Af^ - A^""^ = 0, the eigenvalues in (|B6|l are 
/i'' = A'^ B'^ . The latter vanishes for n — \, while for 
K = one has A" = B'^ . Their contribution can be regrouped 
leading to 



In 1 



|| [i?"(rt/) - B"m 



9k_ 

J^2 



(r"{TU) - i?"(0) 



(122) 



The limit of T — > can then be taken unambiguously, i.e. 
independent of the sign of U: 



1 



(123) 

This gives the final result for the characteristic function of the 
probability distribution of the center of mass u{w = 0) of the 
interface, to lowest order in e = due — d: 



A2 



2m4 
2 



(124) 



In 1 - i?"'(0+)2^A2 + R"'(0+f^X^ 



B. 2-point probabilities and avalanche-size distribution 



1. General considerations 



Following similar arguments as in Section IV A 1 the 2- 
point generating function can be computed as 



g/^ X:c[u(x;w/2)-w~u{x;-w/2)] 
— p/^ Xxlu{x;w/2)-w-u{x;-w/2)]^'^'^° p-ri[u^-*" 



(125) 



In the second line, which is exact, 



denote the s addle - 



point solutions obtained from r[7i. i>]. Here r[u, u] — 
S[u,v] + Ti[u,v] + O(e^) is the effective action of iS[u,?;] 
given by (21 1, setting A = 0. The third line is only correct to 
1-loop order (i.e. lowest order in e = duc — d), and requires the 
evaluation of Fi at the tree-level saddle point; hence we can 



set u 



A 



, A, tree 



, and V 



,,A,trcc 



, A, tree 



Note that this 



symmetry property carries over to the effective action, hence 



one can set from the outset T[u 



A ,,Al 



An im- 



portant property is that the dependence on w of r[w, v] is the 
same as the one of S[u, v], i.e. o nly t hrough the elastic en- 



ergy in (21 1. We can thus derive ( 125 1 w.rt w, to obtain two 



alternative expressions: 



ft pj^ ^3:[n{x;w/2) — w—u{x; — iu/2)] 



(126) 



(127) 



The first one ( |126| l is exact in terms of the exact saddle point, 
hence only explicit derivatives w.r.t. w are needed. The sec- 
ond is expressed in terms of the tree saddle point and is true 
to 1-loop order. Fi has no explicit w dependence, hence the 
last derivative acts only on the dependence on w of the tree 
solution (emphasized in the notation). 



2. Computation o/Fi 

The general expression of Fi [u, v] is given in Appendix [C| 
Again we restrict to uniform A^r = A and uniform Wx ~ w, 
thus we only need the expression for Fi [u, —u] at the uniform 



This result is in agreement with Eq. (G2) in ll20l . At depin- 
ning, the distribution is different, see Eq. (42) of ISOl . 



^ This is because changing w results in a sliift {wl^(x) — Wa{x) + 
w/2, ui^(x) — >■ w^{x) — ui/2) in the generating functional , w^] 

associated to S[u,v]. 
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tree saddle point. Dropping the superscripts A and tree, this is 
Ti[u,-u] 1. 



Trln .g-i^„bll 



1. 
h 



Trln(.g-i<5abll 



Tr 



with 



pK JifK 

^ab ^^^ab 



1 
1 



(128) 



(129) 



(130) 



^ab — 



5ab{-R"{Ual) + R"{Uac) ~ R" {Ua + U^) 



+ R"{Ua+Uc))+KR"{Uab) 
^KR"{Ua + Ub) ■ 



(131) 
(132) 



Here Tr refers to a trace over replica and u, v indices (space), 
and c to c 7^ 1. More explicitly, the matrices and have 
for any n again each four distinct components (denoting a, h 
indices different from 1), 



1 r 

T 



i?"(T[/) - B!'{2u{) + i?"(2Mi - TU) 



+ {k-1)R"{Q) 



/If 



^R"{TU) 



1 



R"{0) - R"{TU) - R"{2ui - TU) 



R"(2ui - 2TU) 



and 



= -i?"(0) 



Pfi = -R!'{2u,) 

Pab = 5abP^ + 



T 



R"{2ui-TU) 



Pc = 



T 



R"{2ui - 2TU) 



(133) 

(134) 
(135) 

(136) 
(137) 

(138) 

(139) 

(140) 
(141) 

(142) 



We recall that U and ui are solutions of ([34] 



The matrix is diagonalized in Appendix [C] for n = 0. 
It has 

(i) with multiplicity —4 the eigenvalue (since = 0). 



Since in ( 136 1 does not depend on k, the contribution of 



these eigenvalues to ( 128 1 cancels between k — 1 and k — 0. 

(ii) 2 eigenvalues each for a — ±1, of the form /i = 5 (Acri 
A^Ba), where p.^, A^, and will be calculated below. 
Regrouping in Eq. ( 128 1 we get 



ri[u, 



E 



a=±l 



In I [9k 
In 



'+lf^ir-\AiBi 



quadratic part 



2 _ ^ a"b" 



(143) 



We find using formula (B15 1, (B16 1 of appendixfB 



i?"(2ui - 2TU) - R"{2ui) (1 - ctk) 



T 

2UR'"{2ui){aK-l) +0{T) 



(144) 



Al = Bl 



Al = {cT- l)U'^R""{2ui)T + 0(r2) (145) 
^i?"(0) + (Ti?"(2ui)l + 0(T°) (146) 



2i?"(2ui - TC/) - R!'{2ui - 2TU) 
+ 2R"{TU) ~ R"{2ui) - 2i?"(0) 
2R"{TU) - 2R"{0) + 0{T'^) 



(147) 



In Eq. ( 143 1 appears the product A^B^, which, contrary to 



each factor, has no ambiguity at T = 0, 



A^Bl = AR"'{Q+fU^ 



Putting everything together we obtain 



(148) 



^ ^ ^ {sin^fc ' + In ([.9fe ' - 2UR!"{2u,)]' - 2U^ R"" {2u^)[R!' {2u,) - i?"(0)] 
- 21n {jg^^ - UR"'{2ui)f - R"'{Q+fU^^ - quadratic partj , 



(149) 
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where subtraction of terms quadratic in U, or equivalently in Rearranging gives 



gfc, is indicated. Inserting into formula (125 1 this gives the 
correction to the tree expression for the 2-point generating 
function for arbitrary w, the distance between the points. It is 
expressed in terms of U and ui which, we recall, are solutions 
of (34 1. Since ui — > 0+ as w ^ 0+ we check that Fi [u, — m] 
indeed vanishes as w — > 0+ as it should from ( 125 i. The lin- 



ear term in w contains the information about avalanches as we 
discuss now. 



3. Limit of w — > 0"*" and generating function of avalanche 
moments 



We now want to use formula ( |127| l, i.e. 



M)-i.O+ 



(150) 



In Eq. \\A9\ the dependence on w is only contained in U and 
Ui. Since ri[it, — u] vanishes for any [/ as w — > 0+, it is 
of the form ri[u, —u] — wf{U) + 0{w^); hence the depen- 
dence of U on w is not needed and we can consider U as its 
w = 0+ H mit, i .e the solution of m^f/ - i?"'(0+)f/^ = A. 
To obtain (150i we thus replace ui — yw/2, expand to 
linear order in w using the w = 0+ limit given in (40i, 



y=[^- ^- We find i = Z - A, with 

Z = Zt,.eo + yi?""(0+)i?"'(0+)x 



(151) 



k i-2C/i?"'(0+)]2 - 2[/i?"'(0+) 



with Zt^cc — m^U. Specifying to 
ing k km yields 



-m^, and rescal- 



Z — Ztrc 



R""{Q+)md-^Zu 

1 — ISfn^trcc 



(fc2 



1 IZi-iQcSm)"^ 
1 



1 



fc2 + 1 — 2Zti-QcSm fc2 + 1 



(152) 



where the terms appear in the same order as in ( 15 1 1. We have 
abbreviated the characteristic scale of avalanches, 

i?"'(0+) 



(153) 



already introduced in (14i. i?""(0+) ^ e is the small expan- 



sion paramet er, an d as indicated subleading terms are of order 



Equation ( 152 1 has the form 



{Z - Zt,.ee)(l - 25„Zt,.eo) - e6Z{Zt,,,) + 0{e^) , (154) 
and since Z — Zucc ^ it can be rewritten as 

(Z-Zt,.ee)[l-5™(Z + Zt,ee)] ^ 6 S Z (Z) + 0(6^) . (155) 



Z = S,nZ^ + Ztrcc(l - S„,Zt,cc) + eSZ{Z) + 0{e^) 



= S,nZ^ + X + eSZ{Z) + 0{e^) 
Explicitly, this is 

Z^X + SrnZ^ + R"'\0+)m'^-^Zx 



(156) 



(fc2 + 1 - 2ZSm)^ 
^ 1 

k'^ -\- \ 'ZZSjji 
Z S„i 



1 



fc2 + l 



(l + fc2)2 



(157) 



We see that ZS,n always appear together. It is therefore useful 
to introduce the dimensionless function Z of the dimension- 
less argument \Sm, 



Z(A5™) ZiX)S,n . 

Inserting into the above equation yields 

Z=X + Z^ + e/2i?""(0+)m'^-S 
1 



(158) 



2 ■J k 



(fc2 



1 - 2Z)2 
Z 



z 



(l + fc2)2 

where the combination 

a := e/2i?""(0+)m'^-4 



k^ + l-2Z k^ + l 
72 

+ 0(62), 



R""iO^ 



(159) 



(160) 



is the fourth derivative of the rescaled renormalized disor- 



der, and I2, defined in equation (67i, is the (dimensionless) 
1-loop integral used to eliminate the normalization of J^,. 
The result ( |159| l is equivalent to Eq. (151) of |20|, noting 
that the force-force correlator used in Eq. (143) of [201 is 
A"(0) = -i?""(0+). 

The generalization to a more general elastic kernel is 
straightforward, and can be obtained replacing k'^ + 1 g^^ 
as detailed in Appendix E of EOl . and for the contact-line ex- 
periment in 145] . 



C. Avalanche-size distribution 

Let us recall the results for the normaUzed probability- 
distribution function p{s), defined in Eq. ( 13 i, as obtained in 



1 19, 20| for standard elasticity; For d > 4, the tree or MP 
result is relevant. It reads 



Pmf(s) 



1 



s/4 



(161) 
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For dimension d smaller than 4, loop corrections are relevant. 
The 1-loop result, obtained by inverse-Laplace transforming 
(|T59l)is 



^exp(c^s-^s'] , (162) 



p{s) 
with exponents 



-^ + ^«-^4(l-C0. (163) 
a 1 

J = l-- = 1 + -(1-Ci)e (164) 

where a = — 1(1 — Ci)^ ™d = 1/3 for the RF class, rel- 
evant to the present study. The constants A, B and C depend 
on e, and must satisfy the normalization conditions 



poo 

/ ds sp(s) = 1 

^0 

/■OO 

/ ds s^p(s) = 2 
■/o 



At first order in e they are 



A = l + -(2~3-fE)a 

o 

B = l-a(l + ^) 
4 

where 7e = 0.577216 ... is Euler's constant. 



(165) 
(166) 

(167) 
(168) 
(169) 



V. UNCORRELATED AVALANCHES: THE BROWNL4N 
FORCE MODEL (BFM) 

A. The BFM model 

In this section we study the Brownian-force model (BFM), 
which corresponds to a Gaussian bare disorder with a force 
correlator in Eq. ^ of 

- R'^iu) = -i?o (0) + a\u\ , R'^\u) = cr sign(w) (170) 

For a point, i.e. in d = 0, V'{u) performs a Brownian motion 
in u. The potential V{u) is thus given by a so-called random 
acceleration process ||5lTl53l . In the present framework we as- 
sume that the distribution of V'{u) has statistical translational 
invariance, hence the model needs a regularization. It can for 
instance be defined in a periodic box V{u + W) = V{u) with 
W OO, the increments V'{ui) — V'{u2) being those of the 
Brownian motion. 



[V'{ui)-V'iu2W ^a\ui-U2\ 



(171) 



The zero mode then has very large fluctuations, i.e. Rq{0) — 
0{W). The generalization to an interface is straightforward 



with V'{u,x) being a set (indexed by x) of mutually un- 
correlated Brownian motions along u. Note that a dynami- 
cal version of this model was studied in the context of non- 
equilibrium depinning [^9' 341. In d = 0, it is known as the 
ABBM model (see |54J for a review). 

A remarkable property of this model defined in the contin- 
uum, is that it appears to be an exact fixed point of the ERG 
in any dimension d, i.e. the renormalized disorder correlator 



R{u) (for its definition see Section IV A 2 and Sections 11 and 
III of Ref. [20 1) remains of the same form as ( |170| l. More 
precisely 

R"'{u) = crsign(M) , R"'{u) = asign(u) , (172) 



where the rescaled disorder correlator was defined in ( |66| l. Its 
flow, i.e. its dependence on to, is given by the ERG equation 

- mdmR!"{u) = (e - C)^"'(^^) + P[R]"'{u) . (173) 

The /^-function, taken for u > 0, contains only higher deriva- 
tives which vanish for ( |172| l, and this to any loop order. This 
property is detailed in Appendix [l] together with a stability 
analysis, which shows that this fixed point is attractive. More 
precisely, it is at least linearly attractive up to 2-loop order 
The roughness exponent for the BFM can be read off from 
Eq. ( |173| l to be C = e = 4 — d. Hence a — Ad<T with 
Ao = 1/4. 

At this stage, this remarkable property is not rigorously es- 
tablished for arbitrary d. In fact, some of the statements have 
to be qualified, see Appendix [l] It should be considered as a 
(quite solid) conjecture. In d = 0, however, there exists some 
theorems, discussed below, which strengthen the case. 

In d = 0, this model has been studied in the context of the 
ID Burgers equation 1361 [37l I4l1 . Let us recall the connec- 
tion. It is convenient to denote space by w and consider the 
time dependent velocity field v = Vt{w) satisfying the Burg- 
ers equation 



dtv - 



(174) 



in the inviscid limit i' — 0^. It is solved via the Cole-Hopf 
transformation 1551 



(175) 



vt{w) ^^[w - u{w)] = V\w) , 

where u{w) realizes the minimum of 
1 



V^(w) — min 



2t 



{u - wY + V{u) 



(176) 



Hence this is exactly the disordered model in d = with the 
(Burgers) "time" t — \/m?, taken at temperature T ^ v /2 — 
0^ (minimization condition), identical to the inviscid limit. 
At initial time t = {),V = V , hence the initial velocity field 
is vt=o{w) — V'{w). The Burgers velocity correlator thus 
equals the renormalized disorder correlator. 



vt{wi)vt{w2) = -R!'{wi - W2) 



(177) 
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with R = Rq att = (i.e. m = oo). 

Until now, these statements were completely general. The 
BFM corresponds to a choice of a random initial velocity 
Vt=o{'w) with the same increments as the Brownia n mo tion. 
As a consequence, R"'{u) is for all times given by ( 172 1. 



B. Shocks in the BFM and Levy processes 

If we admit that there are no loop-corrections for the BFM, 
then we can conjecture that the (improved) tree level (i.e. 
mean field) result is exact for the BFM in any dimension d. 
From the fact that R"'{u) = a (for u > 0), i.e. all higher 
derivatives vanish, the results of Section [HgI then show that, 
for w > Q, 



(178) 



This should hold in any d for the t wo-po int correlation of the 
center-of-mass displacements. In Z {\) takes the tree 

expression Z[\) = Zticc(A) — A from \ Vl\ 



Z{X) 



(179) 



with Sm — fylTn^ = ot^ ■ In Appendix |e] we show, from 
our saddle-point method, that the same holds for an arbitrary 
number of ordered points w\ < W2 < ■ ■ ■ < Wp, 

with Ai = and jii = — ^j- admits a more 

general formulation 



-L-^ J dw ^l' (w)[u(w)-w\ _ L"^ i dw Z(tJi.{w)) 



(181) 



for any function which vanish aiw — ±oo, derived in 
Appendix|F] Inserting = — \i9{w — Wi) one recov- 
ers ([T80]l. 

Let us now make contact with a remarkable set of re- 
sults obtained by Carraro-Duchon and Bertoin for the Burgers 
equation, i.e. the case d = |36, 37|. 

We recall the definition of a (homogeneous) Levy process. 
It is a real random function X{w), continuous on the right 
with a limit on the left, i.e. it can have jumps. It has ho- 
mogeneous and independent increments, i.e. {X{wi^x) — 
X{wiy\i^i^,,,^p are independent random variables for any or- 
dered set wi < W2 < ... < Wp and any p; and for all 
w < w' the law of X{w') — X{w) is the same as the law 
of X{w' — w) — ^(0). Its characteristic function satisfies, for 
w > 0, and cj e iR, 



(182) 



A Levy process is thus fully determined by its Levy exponent 
(j){uj), with 0(0) = 0. More generally. 



for any function u!{w) which vanishes at w = zLoo ^. (182 1 
is recovered using uj{v) = uj9{w — v)6{v). The Levy- 
Khintchine theorem | 56l then establishes that X{w) is a sum 
of a Brownian motion (with drift) and an independent jump 
process, with measure n(s)ds. (We use sans-serif s in or- 
der not do confuse with s ~ S/Sm used earlier.) Here 
we need the case of (i) only positive or zero jumps (resp. 
only negative jumps); (ii) finite first and second moments 
/ max(s, s^)ri(s)ds < cx). In that case 



buj 



1) n{s) ds 



(184) 



5>0 



(The same formula holds with s — > — s for only negative 
jumps). In ( |184| i uj can be taken in a domain of convergence 
which includes Re(cj) > (but usually is larger). 

A remarkable theorem by Carraro and Duchon |36 | estab- 
lishes that if the velocity field Xt{w) — vt{w) of the inviscid 
Burgers equation is a Levy process (with only negative jumps) 
at initial time (with 0'(O) > 0), then (i) it remains a Levy pro- 
cess with only negative jumps for all times; (ii) its associated 
Levy exponent 0t(w) satisfies itself a Burgers equation 



(185) 



We recall in Appendix |F] a simple-minded derivation of this 
formula. Its solution for w > is obtained by inverting 



0t(a; + t(/)o(w)) = (t)Q{uj) 



(186) 



i.e. 0((a;) = (j)o{ht{uj)) where ht{uj) is the inverse function 
of a; — > a; + t(j>{uj). This was applied to the case of the initial 
Brownian velocity 



leading to ll36l[37l 

0t(w) 

n{s) 



00 (w) 



1 + a'^ujt - Vl + 2a'^ujt 



(187) 



2+2 



CO 

1 



\) n(s] ds 



s<0 

1 



ay2^|s|3/- 



-e Sa-^t 



(188) 



(189) 



This is the same law for the shock-size distribution as the 
mean-field result \ \19\ for the interface! 

We can now identify the results from our present method 
with those in d = 0. Since Eq. ( |175| gives vt{w) — [w — 
u{'w)]/t, in d = the process u{wj~~w in the BFM is a Levy 
process with only positive jumps. This is consistent with the 



above, Eq. ( 179 1, noting 



(190) 



-/ dw bj' {w)X{w) _ dw <t>(ui{w)) 



(183) 



' In Bertoin ||37| <f> is called i/j and lu is called q. 
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where we recall t = The result ( 189 1 then gives the 

P{S) of Eqs. ^ with s = -m^S and = 2a. 



To conclude, we conjecture that the BFM model for the in- 
terface in any d has center-of-mass displacements given by a 
Levy process with positive jumps, i.e. perfectly uncorrected 
shocks. In d = this was proven in ll36l[37l . Since we argue 
that for interfaces for more general disorder (i.e. not restricted 
to the BFM model but with shorter-ranged correlations) the 
mean-field theory becomes exact for d > due, we conclude 
that at (and above) the upper critical dimension the BFM be- 
comes a good description (with ^ = 0) and the center of 
mass of the interface undergoes a Levy process. The e expan- 
sion then allows to compute deviations from the independent- 
avalanche properties. 



VI. GENERALIZATION: TREE-LEVEL DIFFERENTIAL 
EQUATION FOR AN ARBITRARY DISORDER 
CORRELATOR 

We have seen in the previous section that in the case where 
I R'" (u) I is a constant, (i) the generating functions for the joint 
probabilities of u{w) — w at an arbitrary number of points 
is easily computed, from the one at 2 points, (ii) the (Levy) 
exponent of the 2-point generating function itself satisfies the 
Burgers equation, as shown by Carraro-Duchon ll36l . 

Here we show an even more striking result: We find a gen- 
eralization of the differential equation, satisfied at tree level, 
to an arbitrary number of points, and for any disorder corre- 
lator R{u). This equation encodes the complete mean-field 
results developed in this paper. Here we show how it arises. 
The question of its solution, and further applications, will be 
examined elsewhere, but for illustration we discuss an explicit 
solution for periodic disorder at the end of this section. 



A. Observable 

The observable of interest is 

„Zt[A] \{w)[u{w)-w] 



(191) 



a slight generaUzation of ( 75 1 since it is now a functional of 



X{w) (hence the notation with a square bracket) and contains 
information about multiple-point correlations. As in Section 
III we work in dimension d = or equivalently (up to the vol- 
ume factor of L'^) we study the center-of-mass displacement 
in any d. We compute it in the (improved) tree approximation, 
i.e. it is the sum of all connected tree graphs (for details of the 



graphical rules see section VIC below): 



Zt[X\ e-''™ '^'^"''["'"'^"'"l 



E 

g 



(192) 



A(to) = — ^ (w) — ^ — - = — ^ — - 



pZtW X(w)lu(w)-w] 



u{w) — w 



Z{X,w) = Z(A,w) - Aw 



TABLE I: Conventions used for the various sources and generating 
functions. 



It is obtained by the expansion of e-^™ Hw)[u(w)-iv]^ where the 
external lines on the top link to the external u{wi) — Wi fields 
at various Wi, following the graphical rules defined in section 
V.C of I.20J . It will be convenient to use the notation, 

n{w)^-X{w)t, t^l/m^ (193) 

and loosely denote by the same symbol 

Zt[17] := Zf[A = . (194) 

It contains information about all n-point cumulants C**^"^ of 
the displacement, or equivalently in d = of the Burgers 
velocity (see table |l] for the conventions used). Zf [57] can be 
expanded as 



oo „ 
n=2 Jwi,...,Wp 



{wp)ci^\wi, ...,Wp) 
(195) 



p=2 

where the cumulants 

{-t)PCj^\wi, ...,Wp)^ [u{wi) - Wi] . . . [u{Wp) - Wpf 

(196) 

were defined in |20|. There we have seen how to calculate 
them at tree level as a sum over all connected tree graphs G as 



in ( 192 1 and obtained them explicitly for p < 4. 



B. Differential equation 

Here we show that Zjil] can be obtained very elegantly 
from a suitable (functional) generalization of the Carraro- 
Duchon equation which naturally sums up all tree graphs in 
the field theory. The idea is to write an evolution equation in 
the variable t; hence we have emphasized the dependence on 
this variable. 

To achieve this, one needs to generalize Zf [Vl] into a func- 
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tional Yf [ri; u] of two variables u and w, defined as 
^ — ' nl 

X / n{wi)...n{Wn)Ct'^\u{wi),...,u{Wn)). 

Hence it depends on a background field u{w), not to be con- 
fused with u{w), the center of mass of the manifold in a given 
disorder realization, the fluctuating field which is averaged 
over. Then the following property holds: 



Yi [n, u] is the solution of the flow equation 

dtYt[n,u] - - 



^—Yt[n,u]-^Yt[n,u] (198) 
ou[w) oil{wj 



with initial condition 
1 



Yt=o[n,u] = ^ I n{w)fl{w')A{u{w)-u{w')) . (199) 



C. Graphical proof 



By definition one has: 



(200) 



By the notation = we mean to use the graphical rules extend- 
ing the ones defined in ||20| Sec. V.C. as follows: 

1. Draw all connected tree diagrams obtained by the ex- 
pansion of e" ■1'^ ^Mt-'[u(w)-w] ^ 

2. Each external point is a contracted variable 
— J^ft{w)t^^[u{w) — w], i.e. will contribute a 
factor of — il.{'w) / 1, hence it does not depend on the 
background field u. 

3. Each dashed line is a disorder correlator, i?(u(wi) — 
u{w2)), with ni derivatives taken w.r.t. u{wi) and 712 
derivatives taken w.r.t. u{w2), where ni and n2 are the 
number of lines entering the left and right vertex respec- 
tively. 

4. Each solid line is a correlation function at zero momen- 
tum, 5q=o — i/ni^ — t. All points connected with 
such a line have the same argument Wi. In the drawing 



of Eq. ( 200 1, we have distinguished external propaga- 



tors, i.e. lines which end in a fl{w) (in green/grey/thin) 
from internal ones (bold, black). The reason is that the 
factor of t on an external line cancels with the factor of 
1/t which comes with each Vl{w). Thus only internal 
lines carry a factor of t. 



5. Once Y[$l, u] has been evaluated, one sets u{w) w 
to get Zt[0], 



u{w)- 



(201) 



In order to allow for a recursion relation, we perform 
this last step only at the end. 

We now show that there exists a recursion relation for deriva- 
tives w.r.t. internal lines: 

Consider 3fY[f2, u]. Since each internal line carries a factor 
of t, graphically this means a sum over all possibilities A4 to 
mark an internal line (here dotted, red), 

dtYt[n,u] 



EE 

M g 



(202) 



One realizes that above and below the marked propagator ap- 
pear functional derivatives of Y[n, u] itself: gJ^^^ Yti'^, u] at 
the top, and — g^^f^^j Yt ['^, u] at the bottom. This implies equa- 
tion ( |198| l. The initial condition ( |199[ ) is then made to recover 
the exact second cumulant C^'^\ i.e. the quadratic term in the 
graphical expansion. 



D. Consequences and particular cases 

The above mean-field differential equation is remarkable 
in several respects. First it allows to c ompu te explicitly the 
n-point function C^"' by integration of ( 198 1 in a small-i ex- 
pansion. One immediately checks that the terms of order fi^ 
and ri** coincide with the expressions (59)-(61) of |20|. 

We now analyze some special cases, for which the above 
equations simplify. Suppose we want to compute the p-point 
expectation 



(203) 



within the tree approximation. We can use the above formal- 
ism with the choice 



n{w) :~ ujiS{w — Wi) 



(204) 



The two functions il{w) and u{w) have been replaced by the 
two sets of discrete variables Wi and Wi. The variation w.rt. 
u(w) gets replaced by the derivative w.r.t. Wi, s.t. one can write 
the recursion relation directly for Zj, 



E 

Q 



I. y ' 



(205) 
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Rule 3]now gives R{wi — W2) instead of R(^u{wi) — u{w2))- 
The mnctional differential Eq. ( |198| l simplifies to an ordinary 
differential equation, 



d 



d 



^ owi duji 

(206) 

The initial condition to this equation is 



1 ^ 



(207) 



Solving Eq. ( 206 1 iteratively in powers of t reproduces again 



the 71-point functions of |20|, Eqs. (59)-(61). Indeed one also 
has, from ( fT95| ) and p04l i: 



tt{{u>i,Wi}) 

n— 2 ii,...,in — l 



(208) 



Formula ( 206 1 simplifies even more for the p = 2 generat- 
ing function expressed as a function of the position-difference. 



gZt(w,t«) ._ ^-flu{w)-u(Q)-w] 

Equation ( 206| l and the initial condition become 

d - d - 

dtZt{^,u) = -—Zt{uj,w)—Zt{uj,w) 
aw ouj 

Zt=o(^,w) -c<.2[A(0)- A(u;)] . 
As an example, we give the solution up to order 

Zt{uj,w) = uj^ A(0) - A{w) 

+2ujhA'{w) A(0) - A(u;) 

5(A(0)- A(w))A'(u;)2 

^2(a(0)- A(u;))^A"(w 



(209) 

(210) 
(211) 

(212) 



This recursion easily reproduces the 6 first connected Kol- 
mogorov cumulants, explicitly calculated in 1201 . Eqs. (62)- 
(66). 

There is an interesting property related to the expansion in 
w: If one writes 



then the equations for the Zp{uj, t), p < n, close, i.e. 

dtzi + zid^zi = , 

dtZ2 + 2z2di^zi + zid^Z2 = , 

dtzz + iz^idujZi + 2z2dujZ2 + zid^z:>, = . 

More generally 

n 

dfZn + ^ qZqdu:Zn_q+i = , 
9=1 



(213) 



(214) 
(215) 
(216) 



(217) 



with initial conditions z„(w, t = 0) 



j2A(")(0+). 



One particular solution of these equations is z„(a;,t) = 
for n > 2. It con-esponds to the BFM (where A(")(0+) = 
for n > 2), discussed in the previous section, and thus to the 



equation ( 1 85 1 originally derived by Carraro and Duchon 1 36 1 . 
It describes a Levy process with exponent zi{uj, t) = (t>t{'^)- 
The initial condition is 

zi{uj,t = 0) = (t>t=o{Lo) = -tj2A'(0+) = . (218) 

zi is uniquely determined by its initial condition, hence we 
can use the result ( |189| l 



, , 1 + 2(TUjt - Vl + 4o-wt 



(219) 



Inserting into (215 1 we find the general solution for Z2, 



Z2 



1 + iuat 



(220) 



This can be seen by introducing a ^ Inw and b = ln(l + 
\/l + Aatuj) in which variables one gets {da + db)z2 = 
j-Z2. The function F{x) is determined by the initial con- 
dition as 



Fix) 



Hence 

Z2{uj,t) 



1 A"(0+) 
'8a?2 A'(0+)2 



2A"(0^ 



(1 + \/l + 4cra;t)2(l + Aujat) 



(221) 



(222) 



Now one can check that zi and Z2 reproduce the terms 0(w) 
and 0{w^) in (64i obtained there by a completely different 
method. For the general case one can determine the z„ recur- 
sively. Their systematic study is left for the future. 



E. Connection with Exact RG equations 

One easily sees that the generalized Carraro-Duchon equa 
tion 



198[ l together with the definition ( |197| l is equivalent to 

(223) 



the following RG equation for the cumulants 

,Wn) 



n 



E ^{p^l)\{q-\)\'' 



p,q,p+q=n+ 



X [C'^'''{wi,W2, ■ . . ,Wp))d^^C'f\wi,Wp+i, . . .,Wn)] ■ 

Here [. . .] means symmetrization over the n variables 
wi , . . . , w„. The summation over p, q is for p,q > 2 in case 
of STS and p,q > 1 in the absence of STS. One can show that 
if C'(^) = at t = 0, it remains so. In that case the equation 
for (7^") involves only C*^"^^'. In appendix [h| we recall the 
Exact RG (ERG) equations in = 0, and show that neglect- 
ing one term in these equations (which corresponds to loop 



corrections) we indeed recover (223 1 which hence appears as 
a tree approximation 
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F. Including loops 



It is shown in Appendix G 2 that Zj [fl] satisfies a more gen- 
eral evolution equation 



dtZt 



51t 



5Zt 



(224) 



This equation is exact (i.e. valid beyond the tree approxima- 
tion) and equivalent to the ERG equations given in Appendix 
[H] Neglecting the first term corresponds to the tree approxi- 



mation, and is equivalent to the recursion of moments (223 1. 



As shown in Appendix |G 2| Eq. \22A) c an also be obtained 
by replacing the tree level equation ( |198 1 by the equivalently 
exact equation 

dtYt [n,u] = - f lim [— ^ T^^t [n,u]] 



w'^w Sfl{'w') 6u{w) 



Su{w 



sn{w) 



Yt[n,u], (225) 



The first term generates all loop corrections. 

We now consider t he B FM. with statistical translation in- 
variance. In appendix |G2| we show that then a solution for Zt 
can be obtained from 



Zf = /t(Woo) + / (l)t{uj{wi),OJoa) 



U!{w) 

dt(l>t{x,y) 
dtftiv) 



dw' ^{w') , Woo = 
(j)t{x,y)di(f)t{x,y) 



(226) 

n{w) (227) 
(228) 
(229) 



where di denotes the partial derivative w.rt. the first argu- 
ment. The initial condition for the BFM is 



(/)t=o(a;,y) 

/t=o(y) 



(TX 
1 



;A(0)y 



axy . 

2 



(230) 
(231) 



The solution of the system ( |228[ )-( [229| with this initial condi- 
tion is 



ftiy) 



1 



l + 2atix 



-Jl + 4at{x- 



a^t^y"^ 



lln(l-tW) + iA„(0)y2 



(232) 



(233) 



The In term corresponds to 1-loop corrections, while for this 
model higher loop contributions identically vanish, as dis- 
cussed in appendices |G 2| and [121 



Here we compute (in d = for simplicity but extension is 
straightforward) the most general 2-point generating function 
using an arbitrary value for a. The calculation is performed 
in Appendix [k| The general result for any function \{w) = 
— on the interval [0^ , 1^] (this is sufficient since u{w) — 
w is periodic) with /i(0) = is 



-Sx 



(234) 



Sx^- I dw[n{w)+m^V{w)]+(TA^ . (235) 
/o 



V{w) 



2<tA 



Act 



2a 



(m2 - 2(7^)2 



A is given by the self-consistent equation 



A^ 



dw 



w?V{w) 
2aA + 2aV{w) 



(236) 



(237) 



From this the 2-point function is obtained by taking /i(w') = 
X9{0 <w' < w): 



^QA[M(tu)-u)-«(0)]^ ^ 



-w 



(238) 
(239) 



V and A both depend on the length of the interval w, and the 
function V{w') introduced above is V{w') — 0(0 < w' < 
w)V . The self-consistent equations are 

Vm^ - 2a AV + aV"^ = -A (240) 
Aim^ - 2a A + 2aV) = wVw? (241) 

to be solved for the branch such that V = + 0(A2) and 
A = + 0(A2). Up to order t^^ , or equivalently A*" or a^ 
this gives, denoting t = l/ra^: 

-Sx^ -{\tfa(w - l)w + 2{\tfa'^t(w - l)w{2w - 1) 
-{\tfah'^(w - l)w{2A{w - l)w + 5) 
+2{XtfaH^{w - l)w(2w - l)(44(w - l)w + 7) 
-2{\tfa^t'^{w-l)w X 

x(52(w - l)w(lA{w - l)w + 5) + 21) 
+0(^l^A^) (242) 



This is in agreement both with o ur previous expansions in ( 64 1 
(up to order w^) and with (212 1. 



VII. CONCLUSION 



G. Periodic case 



In the periodic case in any dimension it is conjectured that 

R"{u) - R"{0) = i?"'(0+)u(l - u). Noting a := R"'{0+), 
we have := i?""(0+) = -2a. 



In this paper we have presented novel and efficient alge- 
braic tools to study multi-point correlations of the displace- 
ment field of an elastic manifold of internal dimension d in 
a random potential, upon variation of an external parameter 
In d = these identify with the correlations of the Burg- 
ers velocity field with random initial conditions (playing the 
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role of the disorder). Such correlations are of interest in the 
field of turbulence. In both cases, they yield the statistics of 
avalanches, i.e. shocks in Burgers. 

The first method uses replica. The saddle-point equations 
obtained at T = resum all tree diagrams and yield among 
others the avalanche-size distribution in the mean-field limit, 
i.e. for d > due- We have then extended this method to com- 
pute the 1-loop corrections. It allowed us to derive the 1- 
loop avalanche-size distribution more systematically than in 
our previous work |20|, providing an independent check of 
the latter. This method has a natural extension to the dy- 
namics, which allows to compute the distribution of veloci- 
ties in an avalanche near the depinning transition [29'l 32ll58l . 
Apart from the avalanche-size distribution, other distributions, 
which were obtained by a resummation of diagrams, as the 
width of an interface 15 9l l60l . or the distribution of critical 
forces at depinning 1301 . should now be obtainable in a purely 
algebraic way. 

The second method arises from the study of the Brownian 
force model (BFM). That model has the unique property that 
"its mean-field treatment is exact", i.e. summation of tree dia- 
gram yields (almost) the exact result. We have argued that this 
property holds in any d. We also proved the stability, i.e. at- 
tractive character of this model under RG to one loop, but we 
believe it to be valid more generally. Ind = this model iden- 
tifies with the Burgers equation with a (stationary) Brownian 
initial condition. We recalled results from the mathematical 
literature: at all times the velocity field remains a Levy pro- 
cess, implying that the shocks are uncorrected. Furthermore 
it was shown that the Levy exponent of this process obeys it- 
self a Burgers evolution equation in time, the Carraro-Duchon 
equation. We then pointed out a more general connection be- 
tween the Carraro-Duchon equation and the mean-field the- 
ory of elastic manifolds, not restricted to the BFM, which al- 
lowed us to: (i) show that avalanches in elastic manifolds at 
and above their upper critical dimension are described by a 
Levy process; (ii) derive a Generalized Carraro-Duchon func- 
tional equation, which is in essence the exact RG equation sat- 
isfied by the mean-field theory, i.e. the sum of all tree graphs. 
This allows in particular to recover very efficiently most of 
the results of the first method at the level of the mean-field 
theory. Extensions including loop corrections were presented, 
but their study was left for the future. They should in principle 
lead to another, maybe more powerful method to study loop 
corrections. 

Both methods presented here have recently been extended 
to a manifold with a A^-component displacement field, and the 
results are presented in |,33i,34J . 



Appendix A: General non-uniform Wx 



In the general case of the 2-point function (45 i at points 
Wx — ±^/(x) we need to analyze Eq. (29 1. Expanding for 



small w, we find 

, 9x'' ivi^') - fi^')] - 2ii""(0+) \yix)\ U{x) . (Al) 
The other equation to be satisfied is 

9^^'U{^') - i?"'(0+) sgn{y{x)) U{xf = A, . (A2) 



These equations can first be studied as an expansion in small 

y{x) = f{x)+2B!"{Q+) f gxx'\fix')\gx'x"Xx" 



+ ... (A3) 

Uix) = 9xx'K' + B"'{Q+)gxx'^gu{f{x')){gx>x'-K"f 

+ . . . . (A4) 

The second-order part of U{x) allows to retrieve the second 
moment of local avalanche sizes, from ( fST) , 



pi{SA) = 2R"\<d+)gxx"9x'x"\f{x")\ 



(A5) 



These equations were studied in |20| in the case f{x) — 1 and 
Aa; — \6{x). In that case y{x) > does not change sign and 



the equation ( A2 1 can be studied separately. More generally 
however, we see from (A3i that if f{x) changes sign, y{x) 



will also change sign (at least for small enough A^,), but not 
necessarily at the same location. A general analysis of these 
equations demands a more thorough study. 



Appendix B: Diagonalization of replica matrices 
1. 1-point formulas 



Consider a replica matrix M as in ( 117i specified by the 



four components Mn, Mia ~ Mai —'■ Mi for a =^ 1, 
Maa Mc + M and Mab =■ M, where a ^ & are two 
arbitrary replica indices distinct from L The eigenspaces can 
be split into two groups: 

(i) a 2-dimensional subspace of vectors of the form 



\ VUJ 



(Bl) 



Acknowledgments 



We are very grateful to V. Vargas for pointing out Ref 
and for a useful discussion on Levy processes. We thank J.L 
Jacobsen for an inspiring discussion on tree resummations. 



There the action of M reduces to a simple 2x2 matrix: 



MV 



Mil {n - l)Mi 
Ml Mr, + {n-l)M 




(B2) 
(B3) 
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(ii) n — 2 eigenvectors associated to the eigenvalue /i = Mc 



with uj 



(p) 



(J 



, = 1, 



(B4) 



, n - 2), the 



{n — l)-vector constracted from the (?i — l)-th root of unity, 

withE;=ic^f -0. 

To summarize, the eigenvalues and multiplicities are for 

n -> 



multiplicity — 2 



(B5) 



fi = - (/i ± V AB) , multiplicity 1 for each sign (B6) 



with 

p. = Mil +Mc- M (B7) 

A = Mil + M - Mc - 2Afi (B8) 

B ^ Mii+M ~-Mc + 2Mi . (B9) 

2. 2-point formulas 

The same analysis can be repeated for the 2n x 2n symmet- 
ric matrix M, with 



ab 



P^b M^b. 



(BIO) 



symmetric 



combination, the 



Since the 2x2 structure is obviously diagonalized by the 

'A . ■ ■ ( 

^ 1 and antisymmetric I ^ 

task reduces to finding the eigenvalues of P'^ ± M'^; hence the 
first —4 eigenvalues are for n — > obtained from (|B5[), 



Mc + Pc 

Mr - Pr 



multiplicity — 2 
multiplicity — 2 



(BID 
(B12) 



The remaining four eigenvalues are according to ( B3 i the 
eigenvalues of the two following 2x2 matrices 



^^^(Mii+Pii -(Mi+Pi) 



Ml + Pi Mc - M + Pc- P 



U- = 



Mil - Pii 
Ml -Pi Mc 



-{Ml 
- M - 



-Pi) 

Pr + P 



These four eigenvalues are 



with cr = ±1, and 

Jla = Mil + Mc - M 
A„ = Mil + M - Mc 
B„ = Mil + M - Mc 



a{Pii +Pc-P) 
2Mi + a{Pii + P - Pc 
■2Mi+a{Pii+P-Pc 



(B13) 
(B14) 

(B15) 



(B16) 

2Pi) 
2Pi) . 



Appendix C: Fi [u, v] 

The general expression of ri[u, v] for 2-point observables 
is an extension of ( |105| i: 



Ti[u,v] = -Trln(.g-M,,ll-Wi,) 



(CI) 



Trln(g-i<5„a-W°,) 



Tr 



1\2 



;0\2i 



ab ab 
ab ab 



(C2) 



with 

W, 



ab^xy rji xy 



Sab'^R"{uacix)) - nE" (Uab{x)) 



Sab ^ R"{Ua{x) - Vc{x)) 



(C3) 



w: 



ab^xy rp 



— m^xy 



5abY,R!'{vac{.x)) - KR"{vab{x)) 



Sab'^P!' {Va{x) - Uc{x)) 



(C4) 



ab^xy 



T 



5.yR!'{ua{x)~v,{x))^WX:^ (C5) 



Appendix D: Diagrammatic representation of 1-Ioop corrections 

Let us recall the graphical interpretation of the (improved) 
tree-level self-consistency equation. As already discussed in 
the text Zmf = A + SmZ^Yp is graphically written as 



MF 




= X + 




(Dl) 



As indicated, the blob denotes Z itself. Note that we work 
in rescaled variables, where A'(0+) — —1, and in order to 
lighten the notation, we count the lower vertex as 1 instead 
of A'(0+) = —1, which explains the change in sign w.rt. 
Eq. 

Let us now consider loop corrections. More details can be 
found in ||20l . A graphical interpretation of the dressed prop- 
agator l/(fc^ + m? — 2SmZ), appearing e.g. in ( 157 1 is 




(D2) 
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The notation on the r.h.s. of the equation is as follows: The 
left vertex of each disorder is at 0, the right one at w. This 
is a graphical representation of the antiferromagnetic rule de- 
scribed in |'20|. The outgoing lines all end in a factor of Z 
which is explicitly drawn. 

There appear two classes of diagrams, corresponding to the 
first and second line (of the integral) in Eq. ( |157| l. The first 
class, denoted Ci in Ii20il . can be written as: 



(D3) 

The diagrams on the second line (of the integral) of Eq. \\59) 
are termed class C2 in |20|. They look like a correction to the 
critical force and can be represented as follows 



1. Discrete calculation 






[A'H-A'(0+)] 



z 



z 



, (k^ + m2 - 2Z) B 



(D4) 



Indeed it should be viewed as a loop of A', of which exactly 
one is expanded in w, leading to a loop with one marked ver- 
tex. This is the vertex drawn above. Note that the double line 
needs at least one A'(w) otherwise it cannot start at and go 



to w as indicated. This leads to the last term in ( D4 1 being 
subtracted. 

If one wants the expression in terms of the renormal- 
ized disorder, one has to subtract the contribution propor- 
tional to Jj, ^fc2^3„i2-)2 , giving an additional term — [A'(it;) — 

2Z 



A'(0+)]/,(^ 



Appendix E: Brownian force model: Many point correlations 

In this Appendix we derive the correlation function of the 
center-of-mass displacement for the BFM model for an arbi- 
trary number of points, thereby giving another derivation of its 
Levy process character discussed in the text. We provide both 
a discrete derivation (of the p-point correlation) and a contin- 
uum one (functional average). To simplify notations, we set 
d = 0, which amounts to omitting the factor of L'', restored 
in the main text. Of course this is achieved within tree-level, 
since we have argued that this be exact for the BFM. 



We start with the tree-level equations p6] ), derived for ar- 
bitrary R{u) and specify them to R"'{u) = (Tsign(u). In this 
section, for notational simplicity we set = 1, R"'{Q^ ) — 
CT = 1 and denote A(0) := -i?"(0). We must solve the fol- 
lowing system of equations for Ui and uic 



Ui + ^sgn(iiii - uij)UiUj = \i . 



Insert the result in 



^T.,{^.-hUi){Ul^-Wi) 



(El) 



(E2) 



We choose the Wi ordered as wi < W2 < ■ ■ ■ < Wn- The 
second equation above implies Ui = A^. Hence we 
can shift uu ~> uu + A(0) J^i ^i^^ eliminate A(0) with- 
out changing the equations. Thus the dependence on A(0) is 
trivial, and we now compute the rest setting A(0) — 0: 



aSZi \i{u{wi)-Wi) 



iA(0)(i:.A.)'e-i5:.,(A.;7,+A,C/,-C/.(7,)|«i.-«u| 



(E3) 



A(0)=0 



The solution of the first equation of (El 1 for n points satisfies 



Ui+l - Ui 



Ki <n-l. 



The second set of equations in (El 1 can be rewritten as 

^ i—l n \ 

A, = {7jl + ^C/, - ^ C/jj , l<i<n. (E4) 

Its solution is 

2[/i = 

2t/, = 
2i7„ = 



^Ai 



VA^ - VAj_i , 2>i>n-\ 

n 

1 + ^ A, - ^/K^, 

3 = i 



(E5) 

(E6) 
(E7) 



with 



A,: = 1+(^A,)V2(^A,- 



Hence we find 



Ui+i - u. 



j=l+l 



Wj+i - Wi 

x/Al 



(E8) 



(E9) 
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This allows to rewrite 



^(A,[/, + XjU^ ~ U,Uj)\uu - uij l (ElO) 

n~\ i n i n 

n— 1 n 

— ESr[(EA.)^ + (i-x/A:) 



j = l j=i+l 
\2 



1=1 V I 



+2(l-yA;)(^A,- ^ A,) 

j = l j=i+l 



Using now twice the definition of it can be rewritten and 
simplified to 



ei:'U^^[nM-m]'""' = e5^(")(^-^-)' (Ell) 



This proves formula ( 180 1 in the main text. 
Consider now Xi — u{wi) — Wi and choose 

Al = -fil + ^t-l, X2 = fil - H2, ■ ■ ■ , A„ = fin-l + ■ 

(E12) 

We then find 



-tree 



ei:r=i' ^^^ix,+^~x,)+^^,{x^+x„) 

i.e. the variables Xi+i — Xi are still independent for fixed pi, 
but are not independent of Xi + X„. However, if one consid- 
ers the rescaled variable {Xi + X„)/-y/A(0), then for large 
A(0) one recovers statistical independence. A similar result 
holds with {X1+X2 + ...+ Xn)/^A{0). 

2. Continuous version 

Let us consider equations ( [60| , ( [6T] l in the main text, and 
specify to the BFM model, with a = i?"'(0+). One must 
solve 

m'^{ui{'w) ~ w) + a / \ui{w) — ui{'w')\U{w') 



A(0) / U{w') 



(E14) 



w?U{w) + a I sigTi{ui{w) - ui(w'))U{w)U{w') 



and insert into 

-5a - 



(E15) 



X{w)- —U{w^[ui{w)-w]. (E16) 



Here and below f = f dw. We now restrict to test func- 
tions such that X{w) = 0, hence U{w) = 0. We define 

U{w) = V'{w) , X{w) = -fi'{w) . (E17) 

Let us assume that fJ.{w) vanishes sufficiently fast at w = 
±00, hence the same holds for V{'w) and no boundary term 
arises in any integration by part. The first equation becomes 

m^[ui{w)-w] (E18) 

+a / sign(ui(w) — ui{w'))ui{w')V{w') — . 

J w' 

Now assume that 

sign(ui(it;) — ui{w')') — sign(w — w') , (E19) 
and take a derivative w.rt. w, leading to 



■m? + 2(tV{w) 



(E20) 



Hence monotonicity holds indeed as long as —1 < -^V{w), 
which we now assume, and discuss below. The second equa- 
tion gives 



V'{w) 



SI 



gn(ui(w) - ui{w'))V' {w') 



= X{w) . 
(E21) 



Using monotonicity ( |E19| l and integration by part yields 

V'{w) [m^ + 2(jV{w)\ = X{w) = -pl'{w) . (E22) 
Hence 

m^V{w) + aV{w)^ = -pi(w) , (E23) 

which can be solved as 

,2 



V{w) 



m 
2^ 



1-4— /iH-1 



(E24) 



Note that m'^V{w) = ~'Z[pi{w)) with Z given in ( 17 1. Inte- 
grating by parts we find 



— S'a — 



2a 



_^'(^) _ —V'{w) 



u{'w) — w] 
V{w) 



= a / V{wY = / Z{pi{w)) . 



rri^ + 2(tV{w) 

(E25) 



Here Z(M) = ^Z(;fj/i) and 



Z{pL) = -{l-2pL-^l- 4/i 



(E26) 



which shows formula ( |181| l in the text. Note that the above 
monotonicity condition for ui{w) is equivalent to < 
1, the usual analyticity domain where the generating function 
is convergent. 
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3. Cumulants of u{w) — w (Burgers velocity) 

It is equivalently interesting to obtain, for the BFM model, 
the expression for the cumulants of the renormalized pinning 
force. These were defined in |20|, Section III B, as 

m^''hih2...hn" = L-("-i)'*(~l)"(7(")(«;i, . . . , u;„) 

(E27) 

(E28) 



hi := u{wi) - Wi 



We will choose wi < . . . < Wn- For simplicity, we use di- 
mensionless units setting m = 1, R'"{0~^) = 1, and d = 0, 
all factors being easily recovered. The lowest moments can 
be computed using the tree-level formula (61) in |20|, which 
should be the exact result for the BFM model according to our 
conjecture, giving the following simple expressions: 

C'''^\wi,w2) ^wi~w2~ R"iO) (E29) 
-C^^^wi,W2,W3) ^ 2 {wi ~ 2w2 + ws) (E30) 
C^'^\wi, W2, w^, W4) — 3! (wi — 3w2 + — 1^4) (E31) 

Higher cumulants have been re-calculated and we find that the 
general result can be written as; 

(-1)"(7(")K, ...,wn) = {n- 1)1 J2(]Zl) (-1)'+'^' 

(E32) 

Equivalently 

i-irCniwi,. . . , Wn) = {n - 1)1 a{l - a)"- 

(E33) 

where the rule a' Wi means to expand in powers of a, 
and to replace the i-th power of a by w^. This formula has 
been checked against the Kolmogorov-cumulants K^{w) := 
{{h2 -hiYY ^ a„(w2 - wi) obtained in (96) of l|20l . 
We have also checked that this result is consistent with the 



result for the n-point generating function (El 1 



Appendix F: Derivation of tlie Carraro-Duchon formula 

In this Appendix we give a physicist's derivation of Eq. 
( |185[ ) entering ( |183[ ). For a mathematical derivation see 1361 . 
We use the notation f = f°° dw, and recall that our con- 
ventions are summarized in table |I] 

Consider the Burgers velocity field to be a Levy process at 
time t. Then the Levy-Khintchine theorem lf56l implies that 



/"OO 

n{w) = —uj'{w) , uj{w) = / dull n{wi) 

J w 

for any function ijj{w) such that aj(±oo) = 0. Below we also 
assume that Q.{w) vanishes (sufficiently fast) at infinity. 

Let us assume that it remains of this form at all times, and 
check that it is correct provided that 0^ satisfies some differ- 
ential equation. To show this, first take dt on both sides and 



use the Burgers equation dtVt{w) + ■^dw[vt{w)'^] = 0. This 
leads to 



(F2) 



J W 



I n'{w)vt{w)^eL'M^')^(^') 



2 

2jJ'^'"hn{w)^' 



We have used that jjJ^a;(w') — JJ^ du;i(5(w— wi) = 



w'); then we obtain 



1 



dtcbt(uj{w)) ^ ^ I / / n\w)0{w-Wi)e{w-W2) 



+- 



X (t>'t{^iwi))(f>'t{uj{w2)) 

n'{w)e{w-wi)(i)t{uj{wi)) 



(F3) 



where we have divided by e-^™' ''^'("('" )). Integration by parts 
leads to 



uj\w)d{w — wi)(/)j(a;(u;i))(/)t {ijj{w)) 



+ \ [ u;'{w)cl^U^M) 



9{w - wi)(t)'t{uj{wi)) — (t)t 



f W J Wi 

1 r 



(t>'t{oj{w)) 

J —00 

<P'ti^{w))(j)t{i^{w)) , 



(F4) 



since uj{w) vanishes at infinity (and (pt vanishes in zero). 
Since this is true for any function a;(w), it implies 



dt(f>t{uj) + (t>tiuj)d^(f)t{uj) = 



(F5) 



This is nothing but Eq. (185 1 
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Appendix G: Beyond the Carraro-Duchon formula: including 

loops 

1. Evolution equation 

We now give the more general evolution equation, valid be- 
yond Levy processes. We define, as in the text. 



(Gl) 



where Zt is a priori an arbitrary functional of r2(u'), and we 
only assume that n{w) vanishes at infinity (i.e. J^^n{w) is 
not necessarily zero). Similar manipulations as above in (|F2[l, 



(F3 I using Burgers' equation yield the exact evolution equa- 



tion 



(G2) 



Inserting (195 1 and expanding in Q, this equation provides 



yet another derivation of the exact RG equations (223 1 for 
the cumulants C'"-*. (See Appendix [h| for a derivation using 
replica). 

To make contact with our equation ( 198| , we note that n{w) 
and u{w) always appear together, and thus 



Su{w) 



u'{w) dw 5Vl{w) 



(G3) 



Setting Y[l^, u] ~^ (possible since we no longer 



derive w.rt. u{w)), inserting this into ( 198 1 and integrating by 



part yields the second term in ( G2 1 



The first term in (G2 1 corresponds to loop corrections and 



to the first term in the more general equation (225 i. One can 
see that these are equivalent as follows: 



lim 
lim 



rYt[r!,u] 



5n{w') 5u{w) 

5 n{w) d S 



(G4) 



rYt[r!,u] 



Replacing u{w) - 
w.r.t. u{w)) yields 



6Cl{w') u'{w) dw 6^{w) 
w (possible since we no longer derive 



lim 



lim fl{w) 



d d \ S 



d 



dw drt{w)^ 



dw dw' J 6il,{w) Sfl{w') 



SQiw) 



(G5) 



Thus ( G4 1 is the 1 -loop correction (in d = 0) to be added to 
(fT98l). 



2. Levy processes and Brownian force model 

We now study some particular solutions of the evolution 
equation (G2i. The first one corresponds to the Levy pro- 
cesses discussed above. Suppose one restricts to the case 
where J n{w) — 0, and where Zt is a function of uj{w) — 

Cdw'niw'), 



(/)t(a;(w)) . 



(G6) 



Using that 



5n(w) 



J dwi 6{w — wi) 



one recovers 



equation (F5 i. Equation ( |G2[ i could thus be used to study de- 
viations from Levy processes. 



For a Levy process, we note that the first term in ( G2 1 van- 
ishes when il{w) = 0, i.e. there are no loop corrections to 
averages of velocity differences, and the tree approximation is 
exact for such observables. 

An interesting generalization, within Levy processes, is to 
allow for il{w) ^ 0, i.e. study observables which involve 
the full velocity and not simply velocity differences. In par- 
ticular, we want to know the full solution for the BFM. In 
the case of discrete p-point correlations of the BFM this was 
done in Appendix |E] A generalization of the Carraro-Duchon 
approach allows to treat that case for continuum observables 
such as Zt [51] . An interesting output is that we will recover 
quite simply the full loop corrections for this model obtained 
via ERG in Appendix [H] 

We define as before Lij(w) = dw' n{w'), but now 
w(— oo) — Woo •= dw' Q{w') may be non-zero, while 



n{w) still vanishes at infinity. We show that ( 178 1 is replaced 
by 

= /t(Woo) + / 0t(w(wi), Woo) • (G7) 

For the integral over wi to be convergent at ±oo we need the 
function (j)t{x, y) to satisfy 



(G8) 



which we assume from now on, and which our solution ( |G20[ ) 
given below satisfies. Upon differentiation this also implies 

920t(O,y) = and di(j)t{y,y) + d2My^y) = 0, which we 
use below^. 



Note that assuming the disorder to be statistically translational invariant 
(STS) can be expressed as a Ward identity 



/ 

•J W 



n'{w) 



5n{w) 



obtained by performing the change of vaiiabl es w 
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(G9) 

-> u) + a in all «) in- 
It implies that if Zt 



tegrals appearin g in Zt . e.g. in its definition i 

is a solution of jG2[ then Zt + F{J^^ ^{w)), where F{y) is an arbi- 
trary function, isalso a solution, with a different initial condition. This 
invariance corresponds to adding the so-called Larkin random force in the 
language of interface pinning. One easily checks that STS is satisfied by 
our ansatz. Inserting JoTJinto \G9\ one finds that it vanishes after integra- 
tion by part because ofjGSl. 
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Let us prove that ( |G7| l is indeed a solution of ( G2 1: Consider now the BFM. The initial condition is 



= /t('^oo) 



(GIO) 



9{w ~ Wi)di(j>t{uj{wi),UJoo) + d2(l>t{uj{wi),UJoo) 



sn{wy 



+ dl(t)t{uj{wi),ujoo) 



(Gil) 



Let us first compute the loop contributions, us- 
ing that = 0, and ^ Woo)] = 
— J7(it;)9iyl(cj(?i;), Wqo)- 



1 

^ ^2 

1 r 



1 r 
2 



ri(w) 5i(/)t(a;(u;),a;oo) + 25i920t('^(w), i^oo) 
9i0t(O,a;oo) + 2a20t(O,a;oo) 

9i0t(O,a;oo) + di(t>t{ Woo, l^oi 



We now compute the tree 



(G12) 
contribution 



A- 



X / ri'(w)0(w - W2)9i0t(a;(w2), Woo) 



ff{uJoc)+ / 92(/>t(w(wi),Woo) 



0t(O, Woo) ~ 0t(w Woo) 



0, 



(G13) 



^ I n'{w)9{w - u;i)6'(w - W2) 



X 9i (?!)t ( wi ) , Woo )<9i (w ( W2 ) , Woo ) 
n{w)d{w - u;i)9i(/)4(a;(iyi), Woo) 

x9i(?!)t(w(w), Woo) 
- / (?!)t(w(w), Woo)9i0t(w(?«),Woo) • (G14) 



Hence we find that if the unknown functions (j)t(x, y) and 
ft{y) satisfy the following two equations, then (G2i is sat- 
isfied: 



dtMy) = ^ di(j)t{0,y) + di(j)t{y,y) 



(G15) 
(G16) 



h=oix,y) = ax - axy 
f,My) = ^A(0)y^ . 



(G17) 
(G18) 



This can be seen by rewriting 

ZoM = J / n{wiMw2)Aiwi - W2) 

= iA(0)( I nf-a I n{wi)^l{w2){wi - W2) 
= Ia{0){ f - a f uj{wi)n{w2) 

^ J Jwi>W2 

+a[uj{wi)n{w2){wi - W2)]2]Z 



] — 00 

-■W2 



(G19) 



Using that uj{oo) = 0, the boundary term is zero and one 
obtains ( pTj i, ^&18[ . 

The solution of the system ( |G15| l-( [G16| l with this initial 
condition is 



<t>t{x,y) 



My) 



2(7t2 



l + 2(7tix 



(G20) 



1 + 4<Tt (a; - I) + c^^t^y^ 



iln(l-tW) + iAo(0)y^ 



(G21) 



One checks that (j)t{x,y) satisfies the conditions (G81, and 
that this recovers the loop corrections obtained by the ERG 
method in (|H8l)- pTT3] l. 



Appendix H: ERG equations and their mean-field version 



It is instructive to compare the equations ( |223| l with the 
known exact RG equations; we restrict here to d = 0. In Ref. 
||25l section IV.A.2 (arXiv version) the ERG equations for the 
cumulant of the renormalized potential S''"''(ti;i, . . . w„) = 

^ c 

( — l)"F(wi) . . . V{wn) were obtained. It was shown that 
the function 

E ^^"n^a„...^.aj (HI) 



n\T 

n ai,...a. 



satisfies 



Hence the S'^"' satisfy at T = 0: 
2at^(")K,...u;„) 

rr>("+l) /' 

+ E 



(H2) 



(H3) 



(p-iy.(q-iy. 

p,q,p+q=n+l ' ' 

X [5'Ja..o(«^l' ■ ■ ■ Wp)S[l^,,,oiwi,Wp+i . . . Wn)] , 
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where we have used that 



in SectionfVll. It reads: 



J2 dlJ'-PHwa„...,Wa,) (H4) 
1 ,...,ap 

Was, • ■ • )'«'ap) 

ai,...,ap 



We now use that C'*^"^(u'i, . . . , Wn) 
i-l)"dwi ■ ■ . 9u,„ S"*^"^ (wi , . . . ,Wn) to obtain 



(H5) 



-I 
2' 



P,<3 



p+q=n+l ' ' 



x[a^j(C'(P)(wi, . . . ,Wp)C'(«'(wi,Wp+i, . . . 



This formula works when the C are continuous functions of 
their arguments'". Hence we see here that if one takes out 
the first term one recovers exactly the mean-field RG equation 



(223 1 of the main text. Hence this provides a further derivation 
of this equation. 

In the STS case the lowest-order ERG equations (including 
loops) are: 



9tC(2)(u;i,W2) = (H6) 
dtC^^\wi,W2,W3) = ~^dy,^C^'^\wi,Wi,W2,W3) 

- ^dni3C^'^\'Wi,W2,W3,W3) 

-d^,C'^^Hwi,W2)C^^\wi,W3) 
-d^,C^^Hw2,W3)C^^\wi,W3) 

-d^,C^^Hwi,W2)C^^\w2,W3) . (H7) 



An exact solution including loop corrections exists for all 
C^") in the STS-Brownian case (stationary BFM discussed 



Here we study only the case with STS (translation invariance). The non- 
STS case, e.g. a two-sided Brownian force landscape starting at zero vio- 
lates this condition at the origin, and requires a special treatment, which is 
left for the future. 



^(1 



K) = 

(wi, W2) = - W2) - a'^t^ 
{wi,W2,W3) = —2ta'^ {wi — 2w2 + ^3) 



— W4) — Qa^t'^ 



UOf^a^iwi - 5w2 



(H8) 
(H9) 
(HIO) 
(Hll) 
(H12) 



■ 6W3 — Aw4 + W5) 
10u>3 — 10w4 + 5w5 — We) 
(H13) 



One has thus a constant part — Cncr^t" with c„ = {n ~ 1)! for 
n even and c„ — for n odd. If one sets c„ to zero one recov- 
ers the expressions in ( |E33| l. The fact that this simple exact 
solution exists in that case is a consequence of the property 
that r = iS in the sense discussed in SectionlTSl 



Appendix I: FRG properties of tlie BFM model 

1. Stability of the BFM fixed point 

Let us express the FRG equation using the rescaled force 
correlator A{u) = —R"{u) defined in (66 1. To one loop (first 
2 lines) and 2 loops (third line) the FRG flow for A'(u), de- 
rived in ifTsHTTll is 



— A'(u) 
dm 



e~C)A'{u) + CuA"{u) (II) 
^3A'{u)A"{u) - A"'{u) \a{u) ~ A(0) 



-9^ 



(A(u)- A(0))(A'(u)2-AA'(0)= 



with A = 1 for the statics (the BFM model studied here) and 
A = —1 for depinning (the ABBM model generalized to an 
interface). In both cases, there is a fixed point corresponding 
to C = e with 



A'(u) = -0- sign(u) , A(0) - A(u) = a\u\ 



(12) 



We note that 1- and 2-loop corrections identically vanish for 
the flow of A' (m) at this fixed point - even for its generaliza- 
tion to a A^-component field 161]. This is in fact more general, 
and we claim it to be true to all loop orders. Indeed, it is easy 
to see that higher loops bring more derivatives, hence vanish- 
ing contributions. For the statics, it can be checked to four 
loops in d = 1251 . Hence we conjecture that this property 
holds for any d. 

Note that some parts of the effective action are flowing. For 
instance, one has, in the BFM model. 



md„A{0) = -eA(0) - a 



^, i.e. A(0) = C - 



(13) 



and 



Hence A(0) Cm 
the constant C = Ao(0) has to be chosen sufficiently large 
This is consistent with the fact that the model, defined with 
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statistical translational invariance, must be defined with a reg- 
ularization, e.g. a periodic box of size much larger than any 
other scale, as discussed in the text. Another regularization 
would be to choose a Brownian force with origin at u = 0, i.e. 
V'{0) = but this leads to a different FRG equation which 
we leave for future investigations. 

Let us now show that the above fixed point is attractive, 
at least for a class of perturbations (defined precisely below) 
which are at most of the same range than the BFM model. It 
thus defines a universality class in any d. The stability analysis 
is performed to first order in e — 4 — d, within the 1-loop FRG 
equation. We look for perturbations of the form 



A'{u) 



(14) 



Physically acceptable solutions for g{u) must vanish at infin- 
ity and have a regular Taylor expansion in powers of |u| at 
u = 0. One then obtains to linear order in g{u) 

-mdrng(u) = {e-Qg{u) + Cug'{u) + 3ag'{u)+g"{u)a\u\ 

(15) 

Using ( = ewe can rescale a — ?> ea, ct — > ea and u — > au. 
Thus one must solve for w > 

— mdmgiu) = (tt + 3)g'{u) + ug" {u) — ~ag{u) . (16) 

As indicated, we search for eigenvalues a, where a > are 
stable and a < unstable modes. We find a general so- 
lution, noting Lf^{x) the generalized Laguerre-i polynomial 
and U (a, 6, z) the confluent hypergeometric function 



g{u) = CiC-"C/(3 - a, 3, u) + C2C-^LI_^{u) 



(17) 



The first solution e~"J7(3 — a, 3, u) behaves as it^^ at small 
u (not physically acceptable), except when a — 3 = 0, 1, ... is 
a positive integer, in which case the two solutions are linearly 
dependent. Hence we set Ci = 0. The remaining solution is 
the second independent function, indexed by a. 



g,{u) ■.^e--Ll_^{u)^L\{-u) 



(18) 



which has a regular Taylor-expansion at it = for all a. It is 
thus a physically acceptable solution. For non-integer a. 



9a{u) 



for u 



oo 



(19) 



Thus for a > this is a long-ranged (attractive) perturbation 
of ( [I2| . The cases of integer a > must be treated sepa- 
rately, because then L?_^{—u) either vanishes identically, or 
is a short-ranged solution (see below). Taking the derivative 
of Eq. ([I6|) w.rt. a with the solution ([I8]l in mind yields 



(110) 



"i9„i [daga{u)] = -a [daga{u)\ - ga{u) 



For a = 1,2, ga{u) vanishes, and daga{u) is a long-range cor- 
related eigenfunction of the RG flow. For a = 3, 4, 5,... we 
can restrict our analysis to the 2-dimensional space spanned 
by ga{u) and daga{u)- It has a Jordan block structure and the 
general solution of the flow equation is 



gain) = ci(m) [dag{u)\ 



C2{m)g{u) 



(111) 



with 

ci(m) 
C2(m) 



TO 

Too 

TO 

Too 



ci(too) 



C2{ma) + In 



Too 



ci(too) . 
^" for u 



(112) 
(113) 



Note that for a = 3,4,5, .. . daga{u) ^ u 
whereas ga {u) is short-ranged, as we discuss now: 



ga=5{u) 



'(3-u) 
'(6-u)(2-u) 



(114) 
(115) 
(116) 



The functions ga{u) for a = 1, 2 vanish. For negative a, there 
are polynomial solutions, consistent with the asymptotic be- 
havior (|I9]|, 



ga=-i(w) = 3 + u 

ffa=-2(u) = ^(2 + u)(6 + u) 



ga=-3{u) 



6 



lOu + 10 



(117) 
(118) 

(119) 



These solutions are unstable and physically unacceptable, 
since they grow stronger than |it| at large u. They correspond 
to models with even longer-ranged correlations than the BFM. 

For instance the leading short-ranged eigenmode a = 3 
reads 

A(0)-A(u) = e|u|+6(l-e-l"l) 

Ml ™ 

hence it has a positive Fourier transform (as long as 6 > — e), 
and thus corresponds to a physical disorder direction. 

The question remains whether we have found the complete 
spectrum for all physically allowed perturbations. We argue 
that this is indeed the case: First of all, we have found a com- 
plete basis for short-ranged perturbations, the functions ga{u) 
for integer a. Functions decaying as a power-law can be ex- 
panded in the basis ga{u) for non-integer a, or daga{u) for 
integer a. As perturbations depending on to, they decay to 0, 
either as a power-law (non-integer a), or with additional loga- 
rithmic corrections (integer a). In conclusion, the BFM fixed 
point is stable w.r.t. perturbations of R"'{u) which decay at 
least as a power law at infinity. 



2. More on loop expansion 

To obtain a deeper understanding of the properties of the 
BFM we must look at its replicated effective action functional 
r[u], with u = {ua{x)}c,^i^ ^ n-.xew^- The statement that the 
(improved) tree level is exact means that for any replica field 

Ua{x), 



T[u] = Sr[u] 



(121) 
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The flow of the exact R{u) was discussed above, with the 
claim that R"'{u) = fTsign(u) does not flow (i.e. is indepen- 
dent of m), while i?"(0) flows, but its flow is unimportant. (It 
is part of the regularization required to define the model). 
Let us examine the meaning and validity of the property 
The exact RG flow (as a function of m in any d) of 



(121 



the p-replica part of r[u], S'''^''("i2,...,ip) was written in 1*251, 
see Eqs. (384), (385) of the arXiv version. We want to study 
the force correlator (corresponding in d = to the Burgers 
velocity) hence look at du^ . . . du^S'^P\ui2,...,ip) = dPS^^l 
By analogy with the second moment, where we found that 
Q3g{2) — /J'" (jQgg jjoi" floyy (away from coinciding argu- 
ments), we want to examine (?^'+^S'*^^\ If we assume that all 
d'^R and higher derivatives vanish, it is clear from these equa- 
tions that the feeding term for this quantity vanishes. The nat- 
ural conclusion is thus that d^^^S^P^ is zero for all p > 3 for 
the BFM model. Hence ( |I21| i holds in the sense of derivatives 
i.e. up to terms 0{uP) in the p-replica terms. These terms 
are called random force terms since they can be set to zero 
by a STS transformation u°-{x) — > u"'{x) + gxx' f{x') where 
f{x) is a (Larkin) random force, coupling linearly to the dis- 
placement field. In c? = this is sufficient to ensure that the 
tree calculation is exact for the BFM model. Recently we also 
showed this property for the dynamics in d = |29|. 

In 0? > one should worry about x (i.e. space) dependent 
fields, i.e. non-local parts of the functionals S'^^\ From Eq. 
( |104[ ), (see also (461) of |25|, arXiv version) where the com- 
plete localH-nonlocal two-replica functional p = 2 is obtained 
from the exact RG equations to order R?, one finds for the 



BFM model 



Ky[Uab] = (T'^glyiiign(uab{x))sign(uab{y)) 



(122) 



If we take a third derivative it vanishes away from the singular 
points. Thus naively there is no nonlocal part for d'^S^'^^ and 
the same will be true for higher p. This leaves open the ques- 
tion of how the derivatives act on the singular points (where 
two replica fields coincide for some values of their argument 
x). We will not attempt to answer this question here, but leave 
it for future research. However, we emphasize that the fact 
that the tree level is exact for "sufficiently reasonable", i.e. 
uniform or nearly uniform, field configurations is sufficient 
for computing e.g. center of mass observables using tree-level 
formulas. This can be seen from (122 1 since if Uab{x) does not 
change sign, i.e. the replica are in partially ordered configura- 
tions, taking a third derivative again gives zero. Hence we 



can safely assume that (121 1 holds in any d for: (i) the needed 
derivatives of the p-replica part; (ii) partially ordered config- 
urations. This is sufficient to argue that tree calculations are 
exact for the BFM model in most applications. 



Appendix J: Toy model: Markovian and Poisson process for 
avalanches 



In this Appendix we describe two simple toy models: (i) 
avalanche positions being a Markov process, (ii) avalanche 
positions and sizes being a totally uncorrected process (Pois- 
son process). 



First consider a Markovian model where the location w„ 
of avalanche n depends only on the previous one, with the 
"waiting time", or interval between avalanche I = w„+i^„ = 
Wn+i — Wn distributed according to a distribution Q{£) with 
d£Q{i) = 1. Given that a first avalanche occurs in wi, the 
probability that the n subsequent ones occur in Y[i=2 i'^i ' ^« + 
dwi] is thus Q{wn,n-i) ■ ■ ■Q{w2.i)dw2 ■ ■ -diVn, also nor- 
malized to unity. Also assume statistical translation invariance 
with a uniform density of avalanches, noted po, hence a given 
avalanche can occur anywhere with the same probability. 

For this model one shows that the probability that the in- 
terval [0, w] with w > contains n avalanches and that their 
positions are < Wi < W2 < . . . < w„ < w, is given for 
ri > 1 by 

PwH^l, ■ ■ ■ ,Wn) ^ J dwodWn+ld{-Wo) (Jl) 



„+i - 



1 A 



Q{wi) Y\_ Q{Wi+lA)Q{w - Wn) , 



i=l 



with Q[w) = d£Q{£). To prove this one notes 
that the probability that the origin belongs to an interval 
of size wi.o is ^^•°^^^^-°'> ^ hence the probability that the 
first positive shock occurs in [wi,wi + dwi], wi > 0, 



is ^ / dwo 



7jvQ{wi). Insert- 



ing the factors 1 — Hi^il^l''^ ^ ^i) + Ojwj — w)] and 



^^'^ expanding one gets (jjl|. The proba 
bilities that there are n avalanches in [0, w] are thus 

p(")= f dwi ...dWnP^i::Hwi,...,Wn) (J2) 



1 



dwi / dwioQ(wio) 



dwi{w — 'Wi)Q{wi) 



(J3) 



These expression are easier written Laplace transformed, and 
one finds for p„(s) — dw e~^^pw^^ 



Puis) 



{l~Q{s)YQ{sY-' forn>l (J4) 



P^{s) = j^{Q{s)-l + s{£)) 



(J5) 



with Q{s) = /g°°dwe-^"'Q(w) and {£) = -Q'(0), 
which satisfy the normalization Yl^=QPn{^) = V'*' i-S- 

The case where avalanches are fully independent events, 
i.e. a Poisson process of density p^), corresponds to Q{w) = 
poe-P^"^ and {£) = l/po- Then Q{s) = and one finds 

Pn{s) = (s+po)"+^ surprisingly 



3'"^ = -Apowre-P"^ 



(J6) 
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for n > 0, i.e. the Poisson distribution for the number of 
shocks in the interval. More precisely one finds 

pt'Kwi, . . .,Wn) = p'^e''""^0{O <Wi<...<Wn<w), 

(J7) 

i.e. a uniform distribution for the positions of the shocks. 

Let us now add information about avalanche sizes. As- 
sume the process contains only positive jumps u{w) — u{0) — 
J2a Sad{0 <Wa < w), then 



Appendix K: Periodic case 



^\[u{w)-u{0)] 



(J8) 



oo „ 

If we further assume that the avalanche sizes are independent 
events uncorrected with their location, 

p^^\wi,...,Wn;Si,...,Sn) 

n 

= p(;)K,...,w;„)nm), (J9) 
where P{S) is a normahzed probability distribution, one finds 

CO 

eM«(-)-«(o)] ^ (e^^^)" (JIO) 

71=0 

(e^^) : J dSP{S)e^^ . (Jll) 



The case of a general Q{w) can be solved in Laplace, 



, 1 [1-Q(.)][(e^^)-1] 



{£), 



1- c 



AS\ 



(J12) 



For the Poisson process, we find that ( J12 1 simplifies into (s 
Po(l - (e^'^)))"\ hence 



^\lu(w)-u{Q)] ^ gPQto/dSP(S)(o^^-l) 

^ «,/dSp(5)(c^^'-l) ^ ,,wZ{\) 



(J 13) 



in agreement with our general result (we have set d — 0). 
Although here we have assumed pq finite (which is the case, 
e.g. in numerical simulations |28 1) the above formula remains 
valid when the total density of shocks is infinite as long as the 
density for a given size p{S) (also noted n{s) in the text) is 
finite and J dS Sp{S) is finite. One then recovers the Levy 
process formula for the case of only positive jumps. For a 
proper mathematical formulation see If36ll37]| . 

We observe that in the Poisson case (for d = 0) Z{X) ^• 
—poL^"^ at large negative A. This is dominated by the prob- 
ability that there is no avalanche in the interval w. This limit 
can also be written as Z{—oo) — —1/(5'). Hence the mean- 
field formula ^VJ\ is valid only for A > — l/S'min, where S'min 
is a typical small-scale cutoff for the avalanche size, as dis- 
cussed in Sec. V.E of Ii20i . 



We now solve the equations ( |60| ) and ( |6T] l for the case where 
R{u) is periodic. It is sufficient to choose X{w) — —p!{w) 
on the interval [0, 1] with ^(0) = Indeed since in the 

periodic case we know that ui (w) — w is periodic of period 1, 
we can write: 



\{w + n) = ~X{w) (K2) 



with dw A(w) — 0. Using ui{w + n) = ui{w) + n one 
finds for w £ [0, 1] that the equations (|60]) and hold where 
all integrals are over w' E [0,1] and U{w) U{w). We 
define again for w € [0, 1] 



\{w) = -p.'{w) , Uiw) = V'{w) . 



(K3) 



Note that since A(0) = A(l) and dwX{w) = one has 
/i(0) — /i(l) and similarly V{0) — V{1). The second equa- 
tion in ([60| gives after integration by part 



V'{w) 



dw'R""{ui{w) - ui{w'))u[iw')V{w') 



(K4) 



The boundary term [R'" {ui{w) — ui{w'))V{w')]q vanishes 
only if the integral goes from 0^ to 1~, i.e. contains the delta 
function at 0, a convention which we use here. Noting that for 
the periodic case 

R"'\u) = i?""(0) + ^ 2i?"'(0+)<5(w - n) (K5) 

n 

we obtain for w £ [0, 1]: 



m2 + i?""(0) / dw'u[{w')V{w') + 2R"'{0+)V{w) 



xV'{w) = -Ai'(w) 
Integration yields 



(K6) 



Viw) m^ + R""{0) / dw'u[{w')V{w') 
+R'"{0+)V{w)^ = -^iiw) . 



(K7) 



(Note that a possible integration constant has been dropped 
since it does not enter any physical observable.) The first 
equation in ([SO]) gives after integration by part 



+ [ dw'R'"{ui{w) ~ui{w'))u[{w')V{w') ^0 . (K8) 

"'0 
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Taking a derivative and using again Eq. ( K5 i, we finally arrive 
at the system of two equations: 



(K9) 



We have defined A = dw' u'{w')Viw'). We can now 
close the system of equations. This gives the result quoted 
in the text with r4 = —2a. The action becomes 



m2 + i?""(0) /g dw'u[{w')V{w') + 2R"'{0+)V{w) 



V{w) 



m2 + i?""(0) / dw'u[{w')V{w') 



+ R'"{0+)V{wf = ~^i{w) 



(KIO) 



Noting a = R"'{ 0) and = R"" (0) we ob tain, multiplying 
u'-^{w) from ( |K9] | with V{w) given by ( |K10| l 



u^{w)V{w) _^ ^ 2aF(u;)][m2 + + aV{w)] ' 

(Kll) 



S'a 



dw 



-^i'{w) ^V'{w) [ui{w)-w] 



/•I 1 

/ dw [fi{w) + m^V{w)] - -TiA^ (K12) 
Jo 2 



as quoted in the text. 
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